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ABSTRACT 



This thesis reports on a simulation study of parametric and nonparametric 
procedures for obtaining confidence intervals for the logarithm of the probability a 
semi-markov process enters a particular state before a fixed time t. Three estimators 
and confidence interval procedures are proposed and compared. The different 
estimators use different amounts of information about the process. The maximum 
likelihood estimator and its normal confidence interval procedure uses the most; the 
estimator based on the empirical distribution function of the observed first passage 
times uses the least. An estimator based on an exponential approximation to the 
survivor function of the first passage time uses an intermediate amount of information; 
confidence intervals for the last estimator are obtained using jackknife and bootstrap 
procedures. The maximum likelihood procedure is the most efficient if the underlying 
model is correct. If the model is not correct the empirical survivor function estimator 
appears to be best for small times and the estimator based on the exponential 
approximation best for large times. 
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I. INTRODUCTION 



A. OBJECTIVES 

Finite state space semi-Markov models Find application in a variety of areas such 
as queueing theory, reliability, and clinical trials. The application oF these models oFten 
centers on the distribution of a first-passage time to a state or a set oF states 
representing for example the lifetime of a system or marking the end of a busy period 
for a server. Suppose that the observations of the path of the semi-Markov process are 
all that is known about the process. 

This thesis reports the results of a simulation experiment to compare various 
parametric and nonparametric estimators of the natural logarithm of the probability a 
semi-Markov process does not enter a particular state before a given Fixed time t. In 
what follows we will use “In probability" for “ natural logarithm of the probability “. 
The speciFic semi- Markov model and estimators considered are given in Chapter II. 
Chapter III contains the details of the simulation experiment and results. Conclusions 
from the study are given in Chapter IV. 



B. SCOPE OF THE THESIS 

The purpose of this thesis is to use simulation to compare estimators and their 
confidence intervals for the In probability a particular semi-Markov process does not 
enter a particular state before a given Fixed time t. The particular semi-Markov model 
and estimators considered are given in Chapter II. A simulation study comparing bias 
and standard errors for these estimators was reported in Gallagher (1986) [Ref. 1]. In 
this thesis, we are primarily interested in comparing confidence interval procedures. 

An estimation procedure which uses the least information about the semi- 
Markov process uses the empirical distribution of the observed First passage times. An 
estimation procedure which uses the most information is to assume a parametric form 
for the sojourn time distribution in each state and a transition matrix to describe the 
transition between states; the parameters of the sojourn time distribution and 
transition probabilities can be estimated using maximum likelihood. An estimation 
procedure requiring less information uses nonparametric estimators of the sojourn time 
distributions and the maximum likelihood estimators for the transition probabilities. 
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In many cases, parametric assumptions concerning the sojourn time distribution 
are difficult to justify. It was demonstrated in Park (1986) [Ref. 2] that incorrect 
parametric assumptions may lead to very biased estimators. Hence, a nonparametric 
estimation procedure may be preferred to a parametric one when actual data is used. 
However, the nonparametric procedure can be expected to be less efficient than a 
parametric one provided the parametric model is correct. 

The thesis is organized as follows. In Chapter II, the nature of the problem is 
described and several parametric and nonparametric confidence interval estimation 
procedure are introduced. In Chapter III, the simulation experiment is described and 
results are given. Conclusions drawn from the simulation study are given in Chapter 
IV. 
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II. NATURE OF THE PROBLEM 



A. DESCRIPTION OF PROBLEM 

The semi-Markov process model used in the simulations to compare the 
estimators is as follows. 

Suppose we observe N individuals. Let X t (i) be the state of the i-th individual 
at time t. We will assume (X t (i), t^ 0} i = 1 , 2 ,...., N, are independent semi- 
Markov processes with three states {0.1,2} having the same distribution as (X. . tS: 
0}. All individuals start at t = 0 in state 1, (X Q = 1). An individual stays in state l for 
a random length of time having distribution function Fj. Upon leaving state 1. the 
process transitions to state 0 with probability 9 and to state 2 with probability 1-9 . 
If the process transitions to state 2. it spends a random length of time there having 
distribution function F,. From state 2, the process transitions to state 1 with 
probability l. State 0 is an absorbing state. Once in it. the individual never leaves. 

For all individuals, the entire path of transitions and sojourn times are observed 
until the time of absorption in state 0 (Fig. 2.1). 

Let 



D = inf (t> 0, X t =0}, 

be the entrance time to state 0 (or the time of death). 

The problem is to estimate the logarithm of the first passage time survivor 
function P (D>t) for fixed time t from the data obtained by observing the N 
individuals. 



B. THEORETICAL DEVELOPMENTS 

Two estimators for P(D>t) will be described in this section. They are the 
maximum likelihood estimator and the asymptotic renewal estimator from a paper by 
P.A. Jacobs. [Ref. 1: p.2-6] 
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Figure 2.1 Transition State Step (Sojourn Time). 



1. Maximum Likelihood Estimate for Continuous Time Markov Chain 

In this subsection, the maximum likelihood estimator for P[D> t) will be given 
for the special case in which the sojourn time in state i is exponentially distributed 
with mean l, p- (i= 1. 2). 

Let R be the number of transitions from 1 to 2 for one individual. The log 
likelihood function for the individual is 

L = R In (1-6) + In 0 + R In p ; (eqn 2.1) 



- ( 1 + R ) In Pj - pjTj - p 2 T, 



where T- (i= 1,2) is the total time spent in state i before death . 
The maximum likelihood estimators are 




l + R 
A 1 + R 




(eqn 2.2) 



(eqn 2.3) 
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A R 

p 2 T 

l 2 


feqn 2.4) 


Further R is geometric {0} so E{R) = (1 - 0)/0. Since the > 


4 individuals are 



independent, the maximum likelihood estimators using the data for all N individuals 
are 



A N 

0 - — 

N + R 


(eqn 2.5) 


a N + R 

Pi " - 

T l 

% 


(eqn 2.6) 


A R 

P 2 = 

T 2 


(eqn 2.7) 



r~> ^ 

where R is the total number of transitions from 1 to 2 for all N individuals and T- is 

1 

the total time spent in state i for the N individuals. For the estimators based on data 
for all N individuals. 



Var {0} = I (0) " 1 


(eqn 2.S) 


Var (p A j} = I (p^ -1 


(eqn 2.9) 


A _ . 

var (p 2 ) = l (p 2 ) 1 

where 


(eqn 2.10) 


1(0) = N E { 0, 


(eqn 2. 1 1) 
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KPP - NE{ 



1 + R 

-7?- ; 



(eqn 2.12) 



I(p 2 ) = NE{ — 2 } (eqn 2.13) 

P2 

with N being the number of individuals. 

Let D be the time of entrance to state 0 for an individual. Fix t and put 

S = S(t) = P {D> t}; 



then 



X.,+ p., X, + p, 0p, 

S(t) = { ■ 2 - 2 EXP (X 2 t) - EXP(Xj t) 1 . - - } At (eqn 2.14) 

. . A, j 

where X,, X., are the roots of the equation 

©P 1 P 2 + y (Pi + P 2 ) + y 2 = 0 . (eqn 2.15) 

A parametric estimator for S(t), the survival function, is 



a ^ /\ A 

P M (D>t}= P: - EXP(X 3 t) - EXP(X t t) ‘ - v ' 



A A AAA 

V-p. , 0 Pi 



x. 



A A 

Xl V^2 



where X t and X., are roots of the equation 



AAA A A , 

0p,p^ + y (p, + p,) + t = o . 



(eqn 2. 16) 



(eqn 2.17) 



Since the maximum likelihood estimators are uncorrelated, the asvmptotic 

A A 

variance of P^ (D> t], which will be denoted by Var^j, is approximately 

Var M ( p M (D > t>|0.p 1 .p 2 } (eqn 2.18) 
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A ^ 

= Var { 9 } (5S c 0 ) " + Var { p j } (cScpj)* + 


Var 


A ^ -> 

^ P 2 } (C'S tPy) 


where 










3S 

30 


5 ^ 3S_ + 3S_ _^2 

0 3 A ^ 30 QAj 30 






(eqn 2.19) 


3S 


3 i 3X 

S c^> 1 ^2 






(eqn 2.20) 


‘ 


3A 1 3p 3A, 3o 1 






3S 


9o ! r 1 .1,: _ 1 


33 


3 6 


~~ 3A, 

- ( ean 2.21) 


1 

CM 

Cl 

<n 

# 


> A 1 _a 2-" a 2- A 1 


5a i 


0 P 2 


0 A «■* 0 P ^ 

£ 2 


33 


-s ®°1 r 3 2 


*- 

s* 


, 0 2 


- l)te Al “] (eqn 2.22) 


Sa i 


V A « ~ A ~ ^ A- " , 2 

12 wJ £ A, 






3S 


S 1 r 2 a ^ u 




ill * 


i ]- e ; 2 - J (eqn 2.23) 


3X 2 * 


A 1 ~ A 2 u 1 " A 2 j a| 




U 2 


3A [ 


- 3 . P, 
1 ^ 






(eqn 2.24) 


30 


p r =2 * 2 h 









dX, 



op. 




'V 2A i 



(eqn 2.25) 



5X 1 _ 0O 1 (eqn 2.26) 

5p 2 ’ P r ? 2 -2A i 
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A 



A 



The asymptotic variance of In {D> t} = Var^f / (P\j (D> t))*. 

2. Asymptotic Renewal Estimator 

This subsection describes an estimation procedure for P (D> t}, which is based 
on an exponential approximation to the survivor function P{D>t}. Details of the 
approximation which is obtained by asymptotic renewal theoretic results is given in 
Jacobs [Ref. 3]. The approximation improves as t becomes larger. 



Let 9 = N / (N + R ) be the maximum likelihood estimator of 9 (equation 

A 

2.5). Let <p*(a) be the empirical transform of the sojourn time in state i ; that is. if 
S,(i), S-,(i) (i) denote all the sojourn times in state i for the N individuals, then 



A 



A M ‘ 

^ (a) = (I'M.,) I EXP (ax S k (i) ). 
k= 1 



(eqn 2.27) 



The asymptotic renewal estimator of the survival distribution [Ref. 3] is 



A AAA 



P A (D>t} = (b/ji) EXP(-k x t). 



(eqn 2.28) 



A 



where K is the solution to the equation 




(eqn 2.29) 




\= lj= 1 



I I (S k (l)+Sj(2)} EXP[i(S k (l)+Sj(2))}. 



(eqn 2.30) 




(eqn 2.31) 



k= l 



and 



AAA A 



b = (9 /k) (p^K). 



(eqn 2.32) 
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Note that when a = 0, the left hand side of the equation (2.29) is - 0 < 0. As 

a increases, the left hand side of the equation increases to co, thus there is a unique 

/» 

solution k. The solution may be found by numerical search. One possible numerical 
search procedure is the golden section search. method. 

3. Empirical Distribution estimate 

Let dj, d.,, , d^- denote the observed times of absorption in state 0 for the N 

individuals. A binomial estimator for the survivor function is 

N 

P B (D>t) = ( 1/N) y l (ti »)(d;) (eqn 2.33) 

k= 1 

where l( t( oo) ( d i> = 1 d i > L 

0 otherwise. 



C. CONFIDENCE INTERVAL PROCEDURE 

Suppose X is a random variable whose probability law depends on an unknown 

parameter 9. Given a random sample of X: Xj, X 2 , x n the two statistics lower (L) 

and upper (U) form a 100(1 - a)% confidence interval for 9 if 

P(L < 9 < U} > 1 - a. 

Procedures to obtain confidence intervals for the point estimates for In P(D>t) 
are given below. 

1. Maximum Likelihood Estimator 

The maximum likelihood estimator is asymptotically normal as the number of 
individuals N becomes large, [Ref. 4]. For a fixed t, the maximum likelihood 
confidence interval for In P[D>t} is 

{ L , U } = lnP M {D>t} (VVar M )/(P M {D<t}). (eqn 2.34) 

A 

where Var^j is obtained by using the maximum likelihood estimator of 9, Pj, and p, in 
equation (2.1 1) - (2.26). Confidence limits that are larger than 0 are set equal to 0. 
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2. Empirical Distribution of First Passage Time ( Binomial and Normal Approx. 
Confidence Interval Methods ) 

a. Binomial Confidence Interval 

Since the N individuals are assumed independent, the estimator Pg[D>t} 

of equation (2.33) has a binomial distribution with N trials and probability of success 

P{D> t}. The IMSL routine BELBIN was used to obtain binomial confidence intervals 

for P(D>t}. For a description of the procedure see [Ref. 5: p.390 - 391]. The binomial 

confidence interval for In P(D>t] is obtained bv taking logarithms of the upper and 

~ A 

lower confidence limits. If the lower confidence limits for Pg [D> t} is equal to 0. it is 
set equal to 0.0001. 

b. Approximate Normal Confidence Interval. 

If the number of individuals N is large, the binomial confidence interval 
procedure can be approximated by a normal confidence interval procedure as. follows, 
[Ref. 6: p.99 - 100 , p.954 - 955]. The interval is 



( L , U } = Pg(D> t} ± z l . a / 2 v / Varg (eqn 2.35) 

where 

Varg = ( 1/N){ Pg[D>t}x( 1 - Pg[D>t} )} 

This interval [ L , U } is an approximate 100(1 - ot)% confidence interval for P[D>t). 
The interval for In P[D>t) is { In L , In U }. If either L or U is less than 0, it is set 
equal to 0.0001. If either L or U is greater than 1, it is set equal to 1. 

3. Asymptotic Renewal Procedure 

In this subsection, two nonparametric methods are described for obtaining 

A 

confidence interval using the estimator P^ (D>t). One is the jackknife. The other is 
the bootstrap. They are described below. 

a. Jackknife Estimation Method 

The jackknife was first introduced by Quenouille (1956) [Ref. 7], for the 
purpose of reducing the estimate bias, and the procedure was later utilized by Tukev 
(195S) [Ref. S], to develop a general method for obtaining approximate confidence 
intervals. 

The basic idea of the jackknife estimation method is to assess the effect of each 
ol the groups into which the data have been divided, not by the results for that 
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group alone, but rather through the effect upon the bodv of data that results 
From omitting that group. The two bases of the jackknife 'are that we make the 
desired calcufation for ail the data, and then, after dividing the data into groups, 
we make the calculations for each of the slightly reduced bodies of data ootained 
by leaving out just one of the groups. [Ref.^j 

The jackknife procedure is as follows: 

A 

(1) Let Y a jj (t) be the estimate P^(D> t} computed using all the data. 

(2) Let Yj (t) be the computed statistic using that data which omits the j^ 1 

subgroup where j = 1, 2, , k. In this thesis, for the case in which the 

number of individuals N equals 20 the j^ 1 subgroup consists of all data 
corresponding to the j^ 1 individual. For the case in which the number of 
individuals N equals 50 the first subgroup consists of all data corresponding to 
the first 5 individuals, the second subgroup of the second 5 individuals etc. 
Some cases for N = 50 were run with the j 1 ^ 1 subgroup consisting of the j^ 1 
individual. The resulting confidence intervals differed little from. those obtained 
by leaving out 5 individuals at time. Computational considerations lead us to 
use the fewer number of subgroups in this case. 

(3) Define the j^ 1 pseudo-value by 

Y*j(t) = k In Y all (t) - (k-l) In Y-(t). (eqn 2.36) 

(4) The jackknife estimator Y*(t) is 

Y*(t)= (1,'k) (Y* (t) +Y* (t) + + Y* (t)}. (eqn 2.37) 

1 2 k 

(5) The jackknife estimator of the variance of Y* (t) is 

k 

S*'(t) = (L(k (k-l))} V (Y,vj(t) - Y,,(t)) 2 . (eqn 2.38) 

j-1 

Tukey(1958) proposes that in a wide variety of problems the k estimated 
pseudo-values can be treated as approximately independent and identically 
distributed random variables [Ref. 8], to obtain the following confidence 
interval procedure. 

(6) The jackknife confidence interval is computed as follows. 
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(eqn 2.39) 



{ L , U } = Y*(t) ± tj. a/2 VS* 2 (t) 

where tp a / 2 is the upper 1- a/2 critical point of the t-distribution with k-1 
' degrees of freedom. The confidence interval given by equation (2.39) is a 
function of the estimated variance. If either confidence limit is greater than 0, 
it is set equal to 0. 
b. Bootstrap Estimation Method 

Efron(1979) introduced the bootstrap method for estimating the 
distribution of a statistic computed from observations, [Ref. 10]. In this thesis, the 
bootstrap was implemented as follows. 

Suppose data are gathered for N individuals. Let Rp(n) be the number of 
transitions from state 1 to state 2 for individual n. Let {Sj,(i)} be the collection of all 
sojourn times in state i for all individuals. A bootstrap replication is generated as 
follows: To generate data for one individual, one observation is 'drawn at random with 
replacement from (R,^(n)}; call the observation r^; rp+- 1 observations are drawn at 
random with replacement from the collection of state 1 sojourn times [Sj.(l)}; Tp 
observations are drawn at random with replacement from the collection of sojourn 
times in state 2, (S^(2)}. This procedure is replicated N times to generate bootstrap 
data for N individuals. The estimator In P,^ {D>t} is computed using the generated 
data. This completes one bootstrap replication; B bootstrap replications are done. 

A 

The B estimates of In P A {D>t}are ordered. A 100x(l-a)% confidence interval is 
constructed using the a/2 and 1 - a/2 quantiles of the bootstrap estimates of In P^ 
{D>t}. If either confidence limit is larger than 0, it is set equal to 0. In the 
simulations the number of bootstrap replications is B = 100. 
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III. ANALYSIS OF THE CONFIDENCE INTERVAL PROCEDURES 



A. DESCRIPTION OF SIMULATION 

A Fortran program is written to simulate the semi-Markov process and compute 
the confidence interval. All simulations are carried out on an IBM 3033AP Computer 
at the Naval Postgraduate School using the LLRANDOM II Random number 
generating package [Ref. 13]. 

The data for the simulated experiments are generated as follows: 

(1) An individual starts in state 1 at time 0. 

(2) A random number with distribution Fj is generated for the sojourn time in 
state 1. 

(3) A uniform random number is generated. 

If it less than 0. the process transitions to state 0 and the data for one individual is 
complete. If it is greater than 9. the process transitions to state 2. A random number 
having distribution function F, is generated. The procedure then returns to step 2. 
Data are generated for N individuals and are collected as follows: 

(1) The passage times to state 0 for each individual; 

(2) The sojourn times in states 1 and 2 for each individual; 

(3) The number of transitions from state 1 to state 2 for each individual. 

For each replication of the simulation the estimates and confidence interval in 
Chapter II are computed for In P(D>t}. Each simulation is replicated 300 times. The 
true In P[D>t} is computed. For each procedure, the number of confidence intervals 
covering the true value is recorded. 

The number of individuals that are too low ( true In P(D>t}>U; the upper 
endpoint of the interval ) and too high ( true In P{D> t} < L; the lower endpoint of the 
interval ) are also recorded. The average length of the confidence interval is computed 
as well as the standard deviation of the lengths. The standard deviation is computed by 
subtracting the mean length from each length, squaring the results, summing over the 
300 lengths and finally divided by 299. 
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B. SIMULATION RESULTS FOR THE EXPONENTIAL MODEL. 



In this section results will be reported for a simulation experiment in which the 
sojourn time distributions in both states are exponential; Fj is exponential with mean 
1/Pj and F 2 is exponential with mean 1/ p,. Some true values for P{D> t} for this case 
can be found in Appendix A. 

For each replication nominal 80%, and 90 % confidence intervals for In P(D> t} 
for various values for t are computed using each procedure of Chapter II. The times 
considered are t =0.5, 1.0, 1.5, 2.0, 3.0, 4.0. 

The confidence interval procedures for In P(D>t} are the binomial confidence 
interval (= BIN ), and it's approximating normal confidence interval (= NOR ), for 
the In fraction of individuals who have not entered state 0 by time t: the maximum 
likelihood confidence interval (= MLE ), jackknife confidence interval for the 
asymptotic renewal estimator (= J.K ), and bootstrap confidence interval for the 
asymptotic renewal estimator ( = ■ B.T ).- For each procedure the number of intervals 
covering the true value of In P (D>t) is recorded as well as the number of intervals 
that are too high or too low. These results are reported in Tables 1. 3, 5 and 7. Next to 
each coverage count are given the corresponding coverage proportion in parenthesis. 
Below the coverage numbers is given a 95% confidence intervals for the coverage rates. 
This interval is computed as follows. If P is the proportion of (l-a)% intervals that 
cover the true value then a 95% confidence interval for the coverage rate is 

{ CL , CU } = P ± 1.96V P (1- P )/300. 

If a (l-a)% confidence interval procedure is performing well, then this interval should 
cover about (l-a)% of the time. In addition, if a 80% {respectively 90%) confidence 
interval procedure is working well, then out of 300 replications between 226 
{respectively 260) and 254 {respectively 280) confidence intervals should cover the true 
value. Simulations for which the numbers of intervals that cover are outside of these 
bounds are given in bold type. In the first column of Tables 1, 3, 5 and 7 the true 
P{D> t} is given in parentheses. 

In Tables 2, 4, 6 and 8 are recorded the average length of the confidence intervals 
for In P{D>t). The estimated standard deviation of the length is below the average 
length in parenthesis. If a procedure is performing well, it should not only have the 
correct coverage rate but also a small average length. The simulation results recorded 
in Tables I - 4 are for a simulation whose parameter values are Pj = 1.0, p,= 10.0 and 
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0 = 0.5. The number of individuals N is set at 20 and 50, representing a low and 
moderate number of individuals. 

Coverage results for N = 50 individuals are presented in Table 1. The binomial 
confidence interval procedure tends to overcover. The maximum likelihood procedure 
has the right coverage for all cases. The two confidence interval procedures using the 
asymptotic renewal estimators undercover for t^ 1.0 and have the right coverage for 
t^ 1.5. When a jackknife confidence interval does not cover it is often because it is too 
high (true value < L ). 

Sample means and variance of the confidence interval lengths, for N = 50 
individuals, are reported in Table 2. The average lengths for the maximum likelihood 
and bootstrap are very’ close for t larger than 1.0. The bootstrap and jackknife 
procedures have larger average confidence interval lengths than those for the maximum 
likelihood estimator but smaller than for the binomial procedure. 

Results for a simulation with N = 20 individuals are given in Tables 3 and 4. 
They are similiar to those in Tables l and 2. However, the average length of the 
intervals is larger in Table 4 than in Table 2 reflecting the smaller number of 
individuals. 

Results of a simulation with the same parameters as the first but with the 9 = 
2/3 are given in Tables 5-8. Results for N = 50 individuals are given in Tables 5 and 
6. The coverage results are in Table 5. The binomial procedure has better coverage 
than in Table 1. The procedure based on the asymptotic renewal estimators have the 
correct coverage only for t > 2.0. The average length of the intervals in Table 6 are in 
general longer than the lengths in Table 2. 

Results for a simulation with 0= 2/3 and N = 20 individuals are given in Table 
7 and 8. The binomial confidence interval tends to overcover when the true probability 
is greater than 0.5241 or less than 0.2753. Once again the confidence intervals based 
on the asymptotic renewal procedure undercover for t < 2.0. The maximum likelihood 
confidence intervals give the correct coverage. 
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TABLE 1 



COVERAGE RATIO (EXP MODEL) 
N = 50, 0^0.5, Pj = 1.0, p,= 10.0 



time C.I (%) 



coverage 



BIN 



NOR 



MLE 



J.K 



B.T 



0.5 80 % 



(.7866) 



90 % 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



17(.06) 

2601.87) 

23(.08) 

(-83,-91) 

51.02) 
288(.96) 
7(.02) 
(.94.. 98) 



38(. 13) 
2391. SO) 
236.08) 
(.75, .84) 

176.06) 
276(.92) 
71.02) 
(.89, ,9o) 



27<,09) 
2456.83) 
28 .09) 
(.77..S6) 

14(.05) 
2741.91) 
12 .04) 
(.88 .95) 



1 196.40) 
1676.56) 

14(.05) 

(.50..61) 

836.2S) 

2096-70) 

8(.03) 

(•64,.75) 



1516-50) 

1491.50) 

01 . 00 ) 
(.44, .55) 

1 10(.37) 
1906.63) 

oi.oov 

(.5S..69) 



1.0 80 % 



(.6203) 



too high 
cover 
too low 
int.val 

too hish 
cover 
too low 
int.val 



28(.09) 
2451.82) 
27(.09) 
(. 17 , . 86 ) 

S6.03) 

282(.94) 

10(.03) 

(.9i;.97) 



2S(.09) 
245 .82) 
27 .09) 
(.77..S6) 

1 26.04) 
2756.92) 
• 13(.04) 
(-S9-.95) 




70(. 23) 
2 166.72) 

146.05 > 
(.67,. 77) 

45(. 15) 
247(.82) 



( 






87) 



806.27) 

2136.71) 

76.02) 
(.66,. 76) 

45(. 15) 
2501.83) 

56.02) 
(.79.. 88) 



1.5 80 % 



(.4891) 



too high 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



196.06) 

262(.87) 

19(.06) 

(.S4..91) 

1 1 ( .04) 
2786.93) 
11 .04) 
(.S9..96) 



143!; 



08) 
2436.8 1) 
34(.ll) 
(.77..S5) 

1 1(.04) 
2786.93) 
11 .04) 
(.90..93) 



276.09) 
245(.82) 
28 .09) 
(.77..S6) 

14(.05) 

2746.91) 



(.88. .95 



536.17) 
235 .78) 
12(.04) 
(.74, .83) 

326. 1 1 
262 ,S7 
66 . 02 ) 
(.84, .91) 



46(. 1 5) 
2336.7S) 

2 1 6 .07 ) 
(.73,.82) 

236.08) 
2666-89) 

11 .04) 
(.S5,.92) 



2.0 80 % 



(.3857) 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



186.06) 

2606.87) 

22(.07) 

(.83,.9l) 

9(.03) 

2771.92) 

14(.05) 

(.S9..95) 



3 1 (. 10 ) 
2336.78) 
36 .12) 
(.73..S2) 

26^°90) 
22 .07) 
(.86.. 93) 



27(.09) 

2456.82) 

2S(.09) 

(.77,.86) 

141.05) 
2741.91) 
1 2( .04) 
(.88,-95) 



476.16) 
240 .80) 
13 .04) 
(.75,.85) 

221.07) 
2701.90) 
86.03) 
(.87, .93) 



316.10) 
243 (.81) 
26 .09) 
(.77,.S5) 

151.05) 
268( .89) 
17 .06) 
(.S6..93) 



3.0 SO % 



(.2399) 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



20(.07) 
2621.87) 
186-06) 
(.S4..9 1) 

13(.04) 
2831.94) 
41.01) 
(.92. .97) 




27(.09) 
2456.82) 
28 .09) 
(.77,.86) 

14(.05) 
275 .92) 
11 .04) 
(.S9..95) 



3 66 . 12 ) 
2501.83) 
14(.05) 
(.79,. 88) 

20(.07) 
274 .91) 
61.02) 
(.88. .95) 



25(.08) 
2401. SO) 
35 .12) 
(.75..S5) 

1 1( .04) 
26 76.89) 
22 .07) 
(.85.93) 



4.0 80 % 



(.1492) 



90 % 



too hieh 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



146.05) 

2781.93) 

81.03) 
(.90, .96) 

61.02) 
2861.95) 
SI .03) 
(.93, .98) 



14(.05) 
253 .84) 
33 .11) 

(.so;.ss) 

6 (. 02 ) 
26l(.87) 
33 .11) 
(.83 .91) 



276 .09) 
245 .82) 
28 .09) 
(.77,.86) 

141.05) 

275(.92) 

111.04) 

(.89.95) 



35 6.12) 
2481.83) 
17(.06) 
(.78,-87) 

191.06) 
274(.9 1) 
71.02) 
(.8S,.9o) 



20(.07) 
2391. SO) 
41 .14) 
(.75,. 84) 

101.03) 
2626.87) 
2S(.09) 
(.84, .91) 
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TABLE 2 



LENGTH OF C.I FOR In P(D> t} (EXPONENTIAL MODEL', 
N = 50, 6 = 0.5, p/= 1.0, p 2 = 10.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.2161 

(0.033) 


0.1902 

(0.032) 


0.0889 

(0.012) 


0.1546 

(0.080) 


0.1267 

(0.038) 




90 % 


0.2709 

(0.042) 


0.2447 

(0.041) 


0.1141 

(0.016) 


0.2048 

(0.106) 


0.1619 

(0.047) 


1.0 


80 % 


0.3188 

(0.049) 


0.2882 

(0.042) 


0.1768 

(0.024) 


0.2022 

(0.066) 


0.1810 

(0.032) 




90 % 


0.4008 

(0.056) 


0.3717 

(0.055) 


0.2271 

(0.032) 


0.2680 

(0.088) 


0.2344 

(0.041) 


1.5 


80 % 


0.4164 

(0.056) 


0.3798 

(0.053) 


0.2648 

(0.037) 


0.2822 

(0.088) 


0.2623 

(0.045) 


- 


90 % 


0.5237 

(0.071) 


0.4919 

(0.079) 


0.3400 

(0.048) 


0.3740 

(0.117) 


0.3398 

(0.061) 


2.0 


80 % 


0.5238 

(0.082) 


0.4802 

(0.063) 


0.3528 

(0.049) 


0.3773 

(0.120) 


0.3522 

(0.065) 




90 % 


0.6584 

(0.100) 


0.6253 

(0.076) 


0.4530 

(0.064) 


0.5001 

(0.159) 


0.4583 

(0.086) 


3.0 


80 % 


0.7496 

(0.131) 


0.6925 

(0.125) 


0.5287 

(0.074) 


0.5S46 

(0.188) 


0.5449 

(0.107) 




90 % 


0.9404 

(0.164) 


0.9173 

(0.178) 


0.6790 

(0.095) 


0.7749 

(0.249) 


0.7059 

(0.140) 


4.0 


80 % 


1.0498 

(0.260) 


0.9973 

(0.311) 


0.7047 

(0.099) 


0.7988 

(0.261) 


0.73S0 

(0.148) 




90 % 


1.3136 

(0.323) 


1.4207 

(0.737) 


0.9049 

(0.127) 


1.0590 

(0.346) 


0.9589 

(0.194) 
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TABLE 3 



COVERAGE RATIO (EXP MODEL) 
X = 20, 9 = 0.5, Pj= 1.0, p 2 = 10.0 



time 



C.I (%; 



coverage 



BIN 



NOR 



MLE 



J.K 



B.T 



0.5 80 % 



(.7866) 



90 % 



too hiah 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



14(.05) 
2721.9 1) 

14(.05) 

(.S7..94) 

4(.01) 

2826.94) 



(.91, .97 



40(. 13) 
2281.76) 

32 .ir 

(•71, .81] 

406.13) 

246(.82) 

1 46 - 05 ) 
(.7S;.S6) 



24(.08) 
251 .84) 
25 .08) 
(.79,-88) 

14(.05) 
2766.92) 
10 .03) 
(-S9..95) 



1066.35) 

174(.58) 

20(.07) 

(.52,64) 

786-26) 

2 1 16.70) 

1 1(.04) 
(.65,-76) 



140(.47) 

1606.53) 

06:00) 

(.48,-59) 

996-33) 

2016.67) 

06-00) 

(.62.-72) 



1.0 80 % 



(.6203) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



17(.06) 
2506.83 
33 .11) 
(.79,-88) 

jm 
1 8( .06 ) 
(-SS..95) 



) 



40(.13 
2276.76 
33(.l 1 . 
(.71,.S1) 

17(.06) 
265(.88 
1 S( .06 ) 
(-S5..92) 



246.08) 

250(.S3) 

26^.09) 

(-79,-88) 

146-05) 

277(.92) 

96.03) 

(.89-95) 



726.24) 
2 1 06.70) 

186.06' 
(.65 .75 



466. 15) 
2456-82) 
91.03) 
(.77.. So) 



666 . 22 ) 

2246.75) 

10(.03) 

(.70,S0) 

336.1 1) 
2596.86) 
86.03) 
(.82. .90) 



1.5 80 % 



(.4891) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



156-05) 

265(.88) 

206-07) 

(-85,-92) 

15(.05) 

278(.93) 

76.02) 

(.8S..9o) 



266.09) 
2266.75) 
48( . 1 6) 
(.70,-80) 

15(.05) 

2656.88 

20(A)! 

(.85,-92 



246.08) 
2506.83) 
26 .09) 
(.79—88) 

14(.05) 

2776.92) 

96.03) 

(.89..90) 



48( . 1 6) 
247 .82) 
56.02) 
(.78. .87) 

33(.l 1) 
266 .89) 
1 (. 00 ) 
(.So. .94) 



40(. 1 3) 
23 66.79) 
24(.08 
(.74..S3) 

1 9( .06) 
2666.89) 

(.S5-'94) 



2.0 80 % 



(.3857) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



1 16-04) 
2676.89) 
22(.07) 
(.85,-93) 

1 1( .04 
280 .93 
96.03) 
(.9 , .96) 



246.08) 
220(.73) 
56(. 19) 
(-68-7S) 

1K.04) 

2676.S9) 

22(.07) 

(.85-93) 




50( . 1 7) 
245(.S3) 

(.7^;°8 2 J) 

28(.09) 

271(.90) 

16.00) 

(.S7„94) 



276.09) 
2406.80) 
33 .11) 
(.75,-85) 

10(.03) 

2676.89) 

23(.08) 

(-S5-93) 



3.0 80 % 



(.2399) 



90 0/ ' 



/O 



too high 
cover 
too low 
int.val 


236-08) 

271(.90) 

66.02)' 

(.87,94) 


231.08' 
2496.83' 
28 .09' 
(.79,87) 


> 23(.08' 

) 2511.84' 

\ 26(.09' 

* (.79,88) 


i 406.13) 

) 252 .84) 

) 86.03) 

(.80,88) 


19(.06) 
23 86.79) 
43 .14) 
(.75,84) 


too high 
cover 
too low 
int.val 


86.03) 

286(.9o) 

61.02) 

(.93,98) 


S6.03) 
264(.8S' 
28 .09' 
(.84-92) 


14(.05 
) 277 .92' 

) 96 .03 ) 

(.89,95) 


) 22(.07) 

) 276 .92) 

2(.01 ) 

- (.89,95) 


91.03) 
260(.87) 
31 .10) 
(.83,91) 



4.0 80 % 



(.1492) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



196.06) 
27 1(.90) 
106.03) 
(.S7,94) 

36.01) 
2876.96) 

66 . 02 ) 
(.93,98) 



191.06) 

2246.75) 

57 .19) 
(.70;.80) 

36.01) 

2406.80) 

2S(.09) 

(-75,-85) 




146.05) 
277 .92) 
91.03) 
(.89,95) 



36(.12) 
253 .84) 
1 16.04) 
(.81,89) 

21(.07) 
277 .92) 
21 . 01 ) 
(.89,95) 



I5(.05) 

2316.77) 

546.18) 

(.72,82) 

71.03) 
26 f ( . 8 / ) 
32 .11) 
(.83,91) 
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TABLE 4 



LENGTH OF C.I FOR In P{D<t} (EXPONENTIAL MODEL) 
N = 20, 0 = 0.5, pj-1.0, p 2 = 10.0 



— 

time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.3696 

(0.092) 


0.3026 

(0.092) 


0.1446 

(0.033) 


0.2083 

(0.079) 


0.1854 

(0.059) 




90 % 


0.4598 

(0.115) 


0.3897 

(0.122) 


0.1858 

(0.042) 


0.2713 

(0.103) 


0.2449 

(0.081) 


1.0 


80 % 


0.5451 

(0.131) 


0.4703 

(0.125) 


0.2877 

(0.066) 


0.3072 

(0.108) 


0.2899 

(0.081) 




90 % 


0.6797 

(0.164) 


0.6132 

(0.168) 


0.3695 

(0.085) 


0.4001 

(0.140) 


0.3816 

(0.110) 


1.5 


80 % 


0.7156 

(0.167) 


0.6320 

(0.159) 


0.4308 

(0.099) 


0.4668 

(0.180) 


0.4361 

(0.130) 




90 % 


0.8928 

(0.208) 


0.8351 

(0.227) 


0.5533 

(0.127) 


0.6079 

(0.235) 


0.5790 

(0.181) 


2.0 


80 % 


0.9150 

(0.273) 


0.8328 

(0.306) 


0.5739 

(0.132) 


0.6452 

(0.258) 


0.5971 

(0.1S6) 




90 % 


1.1412 

(0.340) 


1.1847 

(0.837) 


0.7370 

(0.169) 


0.8403 

(0.336) 


0.7927 

(0.259) 


3.0 


80 % 


1.3762 

(0.567) 


1.3689 

(0.87S) 


0.8601 

(0.198) 


1.0177 

(0.419) 


0.9324 

(0.307) 




90 % 


1.7105 

(0.663) 


i 23?4 

(1-773) 


1.1045 

(0.255) 


1.3253 

(0.546) 


1.2427 

(0.429) 


4.0 


80 % 


2.1097 

(1.197) 


2.4656 

(2.060) 


1.1463 

(0.264) 


1.3966 

(0.584) 


1.2722 

(0.426) 




90 % 


2.5832 

(1.308) 


3.9914 

(2.720) 


1.4720 

(0.339) 


1.8188 

(0.761) 


1.6956 

(0.596) 
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TABLE 5 



COVERAGE RATIO (EXP MODEL) 
N = 50, 6= 2/3, Pj = 1.0, p 2 =10.0 


time C.I (%) 


coverage 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 80 % 

(.7231) 

90 % 


too hish 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 


28(.09) 
249 .83 
23(.08j 
( .79 ,.87) 

13004] 
277092 
10 .03 
(.S9..95) 


50(.17] 
2271.76 
23 .081 
(.71,-8 1 ) 

20007] 
2621.87 
10 .03 
(-84,-9 1 ) 


i 27009] 

i 2481.83 
i 25 .08 

(.78,87) 

1 12004: 

► 273.91 

i 15(.05] 

(.88,94) 


i 173(.58] 

l 12 1(-40) 

1 (A5'.°4 2 6) 

I 1491.50] 

l 14S(.49) 

> 31.01) 

(.44,55) 


1 2461.82) 

i 541.18) 

01.00) 
(.14,22) 

1 2101.70) 

i 901.30) 

01.00) 
(.25,35) 


1.0 80 % 

(.5241) 

90 % 


too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 


20007] 
2511.83 
29 .10^ 
(.79,S8) 

13004' 
279 .93 
8003) 
(.90,96) 


321.111 
239 .80] 
29 .10] 
(.75,84) 

20007] 

271(.90 

91.03) 

(.87-94) 


1 27( .09] 

> 248 .83 

i 25 .08 

(.78,87) 

) 1 11.04] 

i 274.91) 
15(.05, 
(.88,95) 


l 131044) 1501.50) 

I 160(.53) 1461.49) 

l 9003) 4001) 

(.48,59) (.43,54) 

) 961.32) 101(.34) 

l 200(.67) 1981.66) 

» 41.01) 11.00) 

(.61,72) (.61,71) 


1.5 80 % 

(.3799) 

90 % 


too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 


30010) 
239 .80' 
31 .10' 
(.75.. 84) 

10003; 

279 .93 
11 (.04 
(.90, .96] 


30(.io; 

239080' 

31.10' 

(.75,84) 

) 1 5(.05] 

) 267.89 

) 18.06' 

i (.85,93) 


) 271.09) 

) 248083' 

> 25(.08] 

■ (-78..S7) 

1 1 1004' 

) 274.90 

) 15(.05 

i (.88,95) 


> 86029] 

» 206069) 

) 81.03) 

(.63,74) 

) 531.18) 

) 244(.81j 

) 31.01) 

i (.77,86) 


t 88029) 

i 1991.66) 

13 (.04) 
(.61,72) 

1 511.17) 

i 2431.81) 

61.02) 

(.77,85) 


2.0 80 % 

(.2753) 

'90 % 


too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 


19006 
254 .85' 
27 .09“ 
(.81, .89] 

91.03) 
275(.92' 
16 .05' 
(.89, .95] 


) 191.06 

) 237 .79' 

) 44(.15' 

i (.74,84) 

) 265® 

) 26(.09‘ 

i (.85,92; 


) 27(.09' 

) 2481.83' 

) 25(.08* 

i (.78,87) 

1 1 ( . 04 
) 274090 

) 15(.05‘ 

i (.88,95) 


) 611.20] 

) 2281.76' 

) 1 1(-04 

i (.71,.81) 

) 34(.l 0 

) 261 (.87' 

) 5002)' 

i (.83,91) 


1 55(. 18) 

) 2271.76} 

) 18.06) 

i (.71,81) 

) 241.08) 

) 2621.87) 

14 .05) 
i (.84,91) 


3.0 80 % 

(.1446) 

90 % 


too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 


321.11 

252(.S4' 

16(.05‘ 

(.80..88] 

16(.05 
277(.92' 
71.02V 
(.89, .95] 


) 321.10 

) 220(.73; 

) 4S(. 16 

i (.68,78) 

) 61.02) 

) 246(.82] 

48(. 16 
i (.7S,.S6] 


) 27(.09 

) 2481.83' 

) 25(.0S' 

i (.78,87) 

121.04 
) 273.90 

) 15.05' 

i (.88,94) 


) 451.15' 

) 240080 

) 1 5( .05' 

i (.75,85) 

1 2 ^ 
1 


) 25(.08) 

) 2391.80) 

) 36.12) 

i (.75,84) 

) 1 1 ( .04) 

) 264 .88) 

25 .08) 
i (.84,92) 


4.0 SO % 

(.0760) 

90 % 


too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 


35012 
2421. 8L 
23 .08' 
(.76,S5) 

121.04 
2851.95’ 
31.00 
(.93, .97] 


) 351.12' 

1961.65 
69(.23 
) (.60,71) 

) 11.00) 

) 276092' 

23 .08' 
) (.89.95) 


) 27(.09 

) 2481.83' 

) 25(.0S' 

1 (.78,87) 

12004’ 
) 273 .90 

15(.05' 
> (.88,94) 


) 311.10' 

) 251 .84' 

) 18.06' 

i (.79,88) 

) 16005' 

) 2741.91' 

) 101.03' 

i (.88,95) 


) 191.06) 

) 233 .78) 

) 481.16) 

i (-73,82) 

) 8003) 

) 263088) 

) 290 10) 

i (-84.91) 
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TABLE 6 



LENGTH OF C.I FOR In P[D>t] (EXPONENTIAL MODEL) 
N = 50, 0 = 273, p 2 = 1.0, p 2 = 10.0 



time ' 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.2536 

(0.040) 


0.2261 

(0.038) 


0.1196 

(0.017) 


0.2232 

(0.092) 


0.1914 

(0.048) 




90 % 


0.31S4 

(0.050) 


0.2913 

(0.049) 


0.1535 

(0.022) 


0.2958 

(0.122) 


0.2506 

(0.060) 


1.0 


80 % 


0.3848 

(0.056) 


0.3502 

(0.052) 


0.2383 

(0.034) 


0.2618 

(0.082) 


0.2510 

(0.048) 




90 % 


0.4840 

(0.070) 


0.4530 

(0.069) 


0.3060 

(0.044) 


0.3470 

(0.109) 


0.3243 

(0.064) 


1.5 


80 % 


0.5260 

(0.085) 


0.4823 

(0.080) 


0.3570 

(0.051) 


0.3841 

(0.134) 


0.3635 

(0.074) 




90 % 


0.6612 

(0.107) 


0.6282 

(0.107) 


0.4584 

(0.066) 


0.5091 

(0.178) 


0.4783 

(0.102) 


2.0 


80 % 


0.6815 

(0.126) 


0.6284 

(0.120) 


0.4757 

(0.06S) 


0.5387 

(0.196) 


0.4977 

(0.106) 




90 % 


0.8556 

(0.157) 


0.8281 

(0.173) 


0.6108 

(0.088) 


0.7140 

(0.260) 


0.6584 

(0.144) 


3.0 


80 % 


1.0966 

(0.332) 


1.0598 

(0.468) 


0.7131 

(0.102) 


0.8739 

(0.332) 


0.7S67 

(0.175) 




90 % 


1.3718 

(0.414) 


1.5369 

(0.930) 


0.9157 

(0.132) 


1.1583 

(0.440) 


1.0428 

(0.237) 


4.0 


80 % 


1.7374 

(0.773) 


1.9119 

(1.277) 


0.9505 

(0.137) 


1.2188 

(0.476) 


1.0833 

(0.248) 




90 % 


2.1575 

(0.901) 


3.1281 

(1.968) 


1.2205 

(0.176) 


1.6155 

(0.631) 


1.4384 

(0.336) 
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TABLE 7 



COVERAGE RATIO (EXP MODEL) 
N = 20, 6 = 2/3, 1.0, p,= 10.0 



time C.I (%) 



coverage 



BIN 



NOR 



MLE 



J.K 



B.T 



0.5 80 % 

(.7231) 

90 % 



too hiah 
cover 
too low 
int.val 

too hiah 
cover 
too low 
int.val 



21(.07) 

2571.86) 

22(.07) 

(.S2..90) 

31.01) 
282(.94) 
15(.05l 
(.91 .97) 



50(.17) 
228(76) 
22(.07) 
(.71, .SI) 

2 1(.07) 
2641.88) 
15(.05) 
(.84.92) 



33(.l 1) 
2331. 78) 
34(.l 1) 
(-73,-82 ) 

1S(.06) 
268). S9) 
14(.05) 
(.86. .93) 



140(.47) 

155(.52) 

56.02) 
(.46, .57) 



222(.74) 
77< .26) 



l(.OO) 
-2I..3 ' 



I) 



194(. 



36-01) 
(.59, .70) 



1871.62) 

1136.38) 

06.00) 

6-32,-43) 



1.0 80 % 



(.5241) 



90 % 



too hiah 
cover 
too low 
int.val 

too hiah 
cover 
too low 
int.val 



3 1 (. 10 ) 
2311.77) 
381.13) 
(-72,-82) 

14(.05) 
2701.90) 
161.051 
(.87, .93) 



316.10) 
231 .77) 

3S(.13) 

(.70..82) 

321.11) 
253(.84) 

166.05) 

(.80,-88) 




34(.ir 

(.73,-82) 

186.06) 

2686.S9) 

14(.05). 

(.86-93) 



10K.34) 

1936.64) 






) 



696-23) 
2306.77) 
. 11 . 00 ) 
(-72—81) 



1121.37) 
176(.59) 
121.04) 
(.53. .64) 

761.25) 

2171.72) 

71.02) 
(.67.. / 7) 



1.5 SO % 



(.3799) 



90 0/ - 



7 0 



too hiah 
cover 
too low 
int.val 

too hiah 
cover 
too low 
int.val 



321.11) 
2321.77) 
36(. 12) 
(.73,'. 82) 

141.05) 

2736.91) 

13(.04) 

(.SS..94) 



321.11) 

1991.66) 

691.23) 

*(.61,72) 

141.05) 
250(.83) 
36(. 12) 
(-79..SS) 



321.11) 
2341.7S) 
34(. 11) 
(.73,-83) 

1S6.06) 
2681. S9) 
14(.05) 
(.S6..93) 



811.27) 
21 1(.70) 

81.03) 
(.65, .76) 

526. 17) 
2456-82) 

31.01) 
(.77.. 86) 



671.22) 

2071.69) 

26(.09) 

(.64,74) 

366.12) 

241(.80) 

23(.OS) 
(.76.. 85) 



2.0 SO % 



(.2753) 



90 % 



too hiah 
cover 
too low 
int.val 

too hiah 
cover 
too low 
int.val 



216.07) 

2516.84) 

28(.09) 

(.79..8S) 

71.02) 
2871.96) 
61.02 V 
(.93..9S) 



211.07) 
22 1 (.74) 
5S(.19> 
(.69,79) 

71.02) 
265(.S8) 
2S .09) 
(-85.92) 



32(.l 1) 
235(78) 
33(. 1 1) 
(.74..S3) 

1 8(.06) 
2691.90) 
13(.04) 
(.86,-93) 



64(.2l) 

227(76) 

96.03) 

(,7l„8I) 

33(. 1 1) 
262 .87) 
51.02) 
(.83. .91) 



381.13) 
2271.76) 
35(.12) 
(.71. .SI) 

1 5( .05) 
2636.S8) 
22(.07) 
(.S4..92) 



3.0 80 % 



(.1446) 



90 0/ ' 



/O 



too hiah 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



121.04) 
274(.9 1) 
14(.05) 
(.SS,.95) 

41.01) 

282(,94) 

14 .05) 
(.91-97) 



121.04) 
226(74) 
62 .21) 
(.70.. 80) 

4(.on 
234(78) 
62 .21) 
(.73..S3) 



326. 1 1) 
235 .78 
33(.l Ij 
(74.-83) 

1S(.06) 
270 .90) 
12 .04) 
( .87-93) 



501.17) 
2391.80) 
1 1 ( . 04 ) 
(.75,.S4) 

286.09) 
267 .89) 
5(.02) 
(.80..93) 



161.05) 
237(79) 
47(. 16) 
(74,-84) 

91.03) 
263(.5S) 
2S .09) 
(.S4..93) 



4.0 80 % 



(.0760) 



too hiah 
cover 
too low 
int.val 

too hiah 
cover 
too low 
int.val 



21(.07) 
279(-93) 
01.00) 
(.90,. 96) 

41.01) 

2961.98) 

01 . 00 ) 

(.97-99) 



2 16 -07) 
2151.72) 
641.21) 
(.67,77) 

41.01) 
2321.77) 
64(.21 ) 
(.73 .82) 




44(. 15) 
2441. SI) 
121.04) 
(.77,. 86) 

24(.08) 
271 .90) 
56-02V 
(.87,-94) 



15(.05) 
2316.77) 
546. IS) 
(71,-82) 

66 . 02 ) 
2616.87) 
33 .11 
(-S2..92 
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TABLE 8 



LENGTH OF C.I FOR 'n P(D>t} (EXPONENTIAL MODEL) 
N = 20, 0 = 273, p/= 1.0, p 2 = 10.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.4304 

(0.109) 


0.3617 

(0.104) 


0.1981 

(0.053) 


0.3451 

(0.172) 


0.3168 

(0.128) 




90 % 


0.5359 

(0.173) 


0.4680 

(0.139) 


0.2544 

(0.068) 


0.4494 

(0.225) 


0.4302 

(0.183) 


1.0 


SO % 


0.6594 

(0.170) 


0.5788 

(0.161) 


0.3952 

(0.106) 


0.4630 

(0.207) 


0.4394 

(0.164) 




90 % 


0.S227 

(0.212) 


0.7618 

(0.225) 


0.5075 

(0.136) 


0.6030 

(0.298) 


0.5937 

(0.242) 


1.5 


SO % 


0.9426 

(0.323) 


0.8710 

(0.469) 


0.5931 

(0.161) 


0.7228 

(0.399) 


0.6560 

(0.256) 




90 % 


1.1756 

(0.402) 


1.2131 

(0.773) 


0.7616 

(0.207) 


0.9413 

(0.519) 


0.8865 

(0.379) 


2.0 


80 % 


1.2684 

(0.587) 


1.2540 

(0.899) 


0.8300 

(0.222) 


1.0229 

(0.603) 


0.9008 

(0.359) 




90 % 


1.5765 

(0.68S) 


2.0258 

(1.817) 


0.9311 

(0.286) 


1.3322 

(0.785) 


1.2240 

(0.517) 


3.0 


SO % 


2.21SS 

(1.311) 


2.5096 

(2.079) 


1.1951 

(0.400) 


1.6583 

(1.021) 


1.4301 

(0.574) 




90 % 


2.6995 

(1.400) 


4.1815 

(2.807) 


1.5347 

(0.914) 


2.1597 

(1.330) 


1.9457 

(0.830) 


4.0 


SO % 


3.512S 

(1.957) 


3.1967 

(2.724) 


1.6119 

(0.778) 


2.3080 

(1.443) 


1.9704 

(0.796) 




90 % 


4.0702 

(1.893) 


4.8194 

(3.183) 


2.0699 

(0.993) 


3.0057 
( 1.S80) 


2.6794 

(1.140) 



C. ROBUSTNESS 

The robustness of the estimates in Chapter II was studied with respect to an 

incorrect model assumption about the distribution of the sojourn time in state 1. In 

the previous simulations, the maximum likelihood estimator used the known correct 
model. 

The data for the simulation experiment in this section is generated from the 
following three state semi-Markov process: Individuals start in state 1 at t = 0. The 
probability of a transition to state 0 is 9, and to state 2 is 1-0. From state 2. the 

probability of a transition to state l is 1. State 0 is an absorbing state. The sojourn 
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time in state 2 is exponential with mean l/p 2 . The sojcum time in state 1 is the sum of 
two independent exponentials with means l/p A and 1/p,; that is, the sojourn time in 
state 1 has a hypoexponential distribution. The same Fortran program for the 
simulation is used but slightly modified for the above change. The data generated are 
analyzed by the same Fortran programs for each estimator as in the first section. In 
particular, the maximum likelihood estimator assumes the sojourn time in state 1 has 
an exponential distribution rather than the true hypoexponential distribution. 

For the first simulation results reported in Tables 9 to 12, parameter values of p t 
= 2.0, p^ = 10.0, p 3 = 2.0, and 9 = 0.5 are used. Again two different numbers of 
individuals are used; they are 20 and 50. The simulation is replicated 300 times, the 
coverage numbers and the average lengths of the confidence intervals are computed. 
The actual value of the survival function is computed by inverting the Laplace 
transform of the passage time to state 0 for the semi-Markov process. Some computed 
values of the survivor function can be found in Appendix B. 

Table 9 reports coverage results for the case N = 50 individuals. The use of the 
incorrect model for the maximum likelihood estimator leads to the majority of the 
intervals being too low for t = 0.5 and t= 1.0. The binomial confidence interval has a 
slight tendency to overcover. The confidence interval based on the asymptotic renewal 
estimator undercover for t = 0.5 and t= 1.0. For larger values of t they have the 
correct coverage. The results for the jackknife and bootstrap individuals are similar to 
those in Table 1. Table 10 reports the average lengths of the confidence intervals. The 
average lengths for the maximum likelihood procedure are similar to those of Table 2. 
For the other procedures, the average lengths are smaller. 

Results of the simulation experiment for N = 20 individuals are given in Table 
11 and 12. Table 11 shows the coverage results. The maximum likelihood procedure 
undercovers for t = 0.5 and t= 1.0 and overcovers for t= 1.5 and t = 2.0. The binomial 
procedure once again tends to overcover. The procedures based on the asymptotic 
renewal estimator have the correct coverage for all t except t = 0.5. The average lengths 
of the intervals are given in Table 12. 

Tables 13 to 16 report results of simulation of a semi-Markov model with 0= 2/3 
but with the other parameters the same; Pj = 2.0, p,= 10.0, p 3 = 2.0 . 

Results for N = 20 individuals appear in Tables 15 and 16. The confidence 
intervals based on the asymptotic renewal estimator have the correct coverage for all t 
except 0.5 and 1.0. The results are similar to those of Tables 11 and 12. 
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TABLE 9 



COVERAGE RATIO 1HYPOEXP MODEL) 
N = 50, 0 = 0.5, pj = 2.0, p 2 = 10.0, p, = 2.0 



time C.I (%) 



length 



BIN 



NOR 



MLE 



J.K 



B.T 



0.5 80 % 



(.8652) 



90 0/ ' 



7 0 



too hieh 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



181.06) 

2611.87) 

21 (.07) 
(.83,-91) 

1 06.03) 
2831-94) 
7002) 
(.92, .97) 



35(.12) 
244 .81) 
21 .07) 
(.77,S6) 

18(.06) 
275 .92) 

M 




0 (. 00 ) 
61.02) 
294(.98) 
(. 001.04) 



1 30043) 
165(.55) 
51.02) 
(.49, .61) 

99(.33) 
20l)(,67) 
1 ( 100 / 
(.61, .72) 



149(.50) 

150(,50) 

1 ( 200 ) 
(.44. .56) 

1061.35) 
194(.65) 
0 ( 100 ) 
(.59. .70) 



1.0 80 % 



(.6743) 



90 % 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



22C.07) 

260(.87) 

18(.06i 

(.83,-91) 

1 1(.04 
280(.93 
91.03) 
(.9t,.96) 



31.01) 
229(.76 
33 .11 
(.72.81 



22(.07) 
269 .90) 
91.03) 
(.86. .93) 



21.01) 
145(,48) 
153(.51 ) 
(.43,54) 

01.00) 
2 13(,7 1) 

87(,29) 

(.66,76) 



42(. 14) 
2421. 81) 
16(.05) 
(.76,85) 

28(.09) 
263 .88) 
9(.03V 
(.84..9I) 



396.13) 
240 .SO) 
21 (.07) 
(.75,85) 

21(.07) 

2651.88) 

141.05) 

(.85,92) 



1.5 SO % 



(.5157) 



90 0/ - 



7 0 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



271.09) 

254(.85) 

19(.06) 

(.81,89) 

1 1 ( .04) 
276( -92) 
13(.04) 
(.89,95) 



271.09) 

2351.78) 

3S(.13) 

(.74,83) 

216.07) 

2611.87) 

18(.06) 

(.S3..91) 



61.02) 

2421.81) 

52(.17) 

(.76,85) 

41.01) 
2766.93) 
20 .07) 
( . S9..95 ) 



311.10) 
2471.82) 
22 .07 
(.78,87) 

216.07) 
271 .90 
81.02) 
(.87,94) 



251.08) 
2431.81) 
32(.l 1) 
(.77,85) 

196.06) 

2621.87) 

19(.06) 

(.S4..91) 



2.0 80 % 



(.3930) 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



25(.08) 
2531.84) 
22). 07) 
(.S0,S8) 

161.05) 

2741.91) 

10(.03 

(.8S..95) 



256-08) 

2406.S0) 

35(.12) 

(.75,85) 

161.05) 

2621.87) 

22(.07) 

(.S4.91) 



1 36.04) 
2641.88) 
23(.OS) 
(.84.92) 



71.02) 
282(,94 
1 1(.04 
(-91,-97 



291.10) 

250(.S3) 

211.07) 
(.79,S8) 

206.07) 
2691.90) 

1 1(.04) 
(.86,93) 




3.0 80 % 



(.2278) 



90 % 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



25(.08) 

2471.82) 

2S(.09) 

(.78,87) 

17(.06) 

2731.91) 

10(.03) 

(.88,94) 



256.08) 

2201.73) 

55(.18) 

(.68,78) 

176.06) 

2551.85) 

2S(.09) 

(.81..S9) 



341.11) 

2531.84) 

13(.04) 

(.80,88) 

20(.07) 
275 .92 





4.0 80 % 



(.1320) 



/O 



too hieh 
cover 
too low 
int.val 

too hieh 
cover 
too low 
int.val 



11 (.04) 
266(.89) 

231.08) 

(.85,92) 

61.02) 

2851.95) 

91.03) 

(.93,97) 



1 1 (.04) 
2351.78) 
541.18) 
(.74,83) 

61.02) 

2711.90 

23(.08 

(.87,94 



47( . 1 6) 
246 .82) 
71.02) 
(.78,86) 

291.10) 

267(.89) 

41.01) 

(.85,93) 



316.10) 

2441.81) 

251.08) 

(.77,86) 

196.06) 

2701.90) 

111.04) 

(.87.93) 



20(.07) 

23S1.79) 

42(.14) 

(.75,84) 

106.03) 
264 .88) 
261.09 
(.84,92 
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TABLE 10 



LENGTH OF C.I FOR In P'D> t}(HYPOEXP MODEL) 
N = 50,0 - 0.5, p t = l0,p 2 - 10.0, p 3 = 2.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.1688 

(0.029) 


0.1442 

(0.029) 


0.0886 

(0.011) 


0.0913 

(0.042) 


0.0806 

(0.022) 




90 % 


0.2110 

(0.037) 


0.1852 

(0.037) 


0.1138 

(0.014) 


0.1210 

(0.056) 


0.1043 

(0.02S) 


1.0 


SO % 


0.2S43 

(0.043) 


0.2553 

(0.041) 


0.1764 

(0.022) 


0.1666 

(0.053) 


0.1534 

(0.028) 




90 % 


0.3572 

(0.054) 


0.3291 

(0.053) 


0.2265 

(0.028) 


0.2208 

(0.071) 


0.1988 

(0.037) 


1.5 


SO % 


0.3958 

(0.058) 


0.3606 

(0.054) 


0.2641 

(0.033) 


0.2692 

(0.083) 


0.2496 

(0.047) 




90 % 


0.4978 

(0.073) 


0.4666 

(0.072) 


0.3392 

(0.042) 


0.3568 

(0.110) 


0.3240 

(0.064) 


2.0 


SO % 


0.5139 

(0.080) 


0.4710 

(0.075) 


0.3519 

(0.044) 


0.3777 

(0.117) 


0.3506 

(0.067) 




90 % 


0.6461 

(0.100) 


0.6130 

(0.100) 


0.4519 

(0.056) 


0.5006 

(0.155) 


0.4561 

(0.093) 


3.0 


SO % 


0.7S40 

(0.160) 


0.7261 

(0.157) 


0.5274 

(0.066) 


0.5998 

(0.188) 


0.5576 

(0.111) 




90 % 


0.9S32 

(0.199) 


0.9684 

(0.247) 


0.6773 

(0.085) 


0.7950 

(0.249) 


0.7261 

(0.153) 


4.0 


SO % 


1.1592 

(0.334) 


1.1276 

(0.480) 


0.7030 

(1.088) 


0.8242 

(0.260) 


0.7650 

(0.154) 




90 % 


1.4496 

(0.415) 


1.6667 

(1.017) 


0.9027 

(1.113) 


1.0924 

(0.345) 


0.9980 

(0.214) 
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TABLE 11 



COVERAGE RATIO (HYPOEXP MODEL) 
N = 20,6 = 0.5, Pj-2.0, p 2 = 10.0, Pj-2.0 



time 



C.I (%) 



length 



BIN 



NOR 



MLE 



J.K 



B.T 



0.5 SO % 



(.8652). 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



12(.04) 
278 .93) 
10(.03) 
(-90,-96) 

0(.00) 

290(.97) 

10(.03) 

(.95,99) 



59(.20) 
231 .77) 
10 .03) 
(.72..S2) 

37(. 12) 
260(.87) 
31.01) 
(.83. .92) 



01 . 00 ) 

53(.18) 

2471.S2) 

(.13..22) 

0(.00) 
109(.36) 
1911.64) 
(.31 .42 



1051.35 
183(-61 
121.04 
(.55 .67) 



22ofi73j 



1 351.45) 
162(.54) 

31.01) 

(.48..60) 

82(.37) 

2171.72) 

11 . 00 ) 

(-67,-77) 



1.0 80 % 



(.6743) 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



161.05) 
2541.85' 
30(. 10 
(.80,88' 

51.02) 

285(.95 

10(.03 

(.93,97 



44(. 15) 
226 .751 
30 .10) 
(.70,. 80) 

161.05) 
2741.91 
10(.03) 
(.88.-95 > 



1(.00) 
2211 - ' 



78(. 26 ) 
(.69,.79) 

0(.00) 
2661.89 
34(. 1 1 
(.85. .92 



451.15 
2401. SO 
15 .05 
(.75,.85) 

29(. 10) 
263(.88) 
81.02) 
(.84. .91) 



32(. 11) 
2411. SO) 
27(.09) 
(.76,85) 

1 1 (.04) 
271(.90 
18 .06). 
(.S7..94) 



1.5 SO % 



(.5157) 



/O 



too high 
cover 
too low 
int.val 

too high 
cover 
too lOW 
int.val 



211.07) 

26i)(.87 

19(.06 

(-83,-91' 

51.02) 



19(.06) 

(.S9..95) 



2l(.07) 
242 .81) 
37 .12) 
(-76,-85) 

21 (.07) 
2601. 87) 
19(.06) 
(.83, .91) 



SI. 03) 
257(.86) 
35(. 12) 
(.82, .90) 

21 . 01 ) 
283(.94 
151.05 
(.92 .97 



401.13) 

2441.81 

16(.05) 

(-77..S6) 

24(.08) 
273 .91) 
31.01) 
(.88, .94) 



231.08) 
232 .77 
45 .15 
(.73,. 82) 

2&% 
29(. 10) 
(. 841.92) 



2.0 80 % 



(.3930) 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



11 (.04) 
270(.9O 
19(.06 
(-S7..93 

111.04) 
2841.95) 
51.02) 
(.92, .97) 




1 7( .06) 
2601.87) 
23(.OS) 
(-83,-91) 

71.02) 
284(.95) 

91.03) 
1.92..97) 



391.13) 

246(.82) 

15(.05 

(.78,.86) 

19(.06) 

277(.92) 

4i:oij 

(.S9,.9o) 



18(.06) 

232(.77 

50(.17) 

(.73,.82) 

71.02) 
263(.S8 
30(. 10 
(.84 .91 



3.0 80 % 



(.2278) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



(.04) 


1 3( 


.04) 


311 


i .91) 


2491 


’.83) 


254( 


1 .05) 


38 


M3) 


15 


:.94) 


(.79, 


;.87) 


(.81 . 



51.02) 
28 1 (.94) 

141 . 05 ) 
(.9 ! ..96) 



51.02) 

257(.8o) 

38(.13) 
(.82 .90) 



171.06) 
279(.93) 
4i:on 
(.90,. 96) 



37(. 12) 
2511.84) 
12 .04) 
(.79,.88) 

19(.06) 
277 .92 
41.01) 
(.S9..95) 



171.06) 
2301.77 
53 .18) 
(-72,-81) 

61.02) 
26l(.8 /) 
33(. 1 1) 
1.83 .91) 



4.0 80 % 



(.1320) 



90 % 



too high 
cover 
too low 
int.vai 


41.01) 

279(.93) 

17(.06) 
(.87, .94) 


4 .01) 
226(.7o) 
70 .23) 
(.70, .80) 


44(.15) 
245 .82) 
11 .04) 
(.79, .88) 


381.13) 
2511.84 
11 .04) 
(.81,89) 


too high 
cover 
too low 
int.val 


41.01) 
2961.98) 
01.00) 
(.97, .99) 


21.01) 
2281.76) 
70 .23) 
.71, .81) 


231.08) 
2731. c » 1) 
41.01) 
(.88, .94) 


191.06) 

277(.92 

41.01) 

(.89,95) 



1 7( .06) 
229 .76) 
54 .18) 
(.72, .82) 

26 ()i°s': ) ) 

(.83^91) 
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TABLE 12 




time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.2891 

(0.077) 


0.2217 

(0.083) 


0.1438 

(0.030) 


0.1405 

(0.056) 


0.1296 

(0.046) 




90 % 


0.3591 

(0.095) 


0.2810 

(0.111) 


0.1846 

(0.038) 


0.1829 

(0.075) 


0.1750 

(0.066) 


1.0 


80 % 


0.4845 

(0.113) 


0.4131 

(0.107) 


0.2861 

(0.059) 


0.2694 

(0.105) 


0.2527 

(0.073) 




90 % 


0.6038 

(0.142) 


0.5365 

(0.143) 


0.3674 

(0.076) 


0.3509 

(0.136) 


0.3314 

(0.102) 


1.5 


80 % 


0.6810 

(0.175) 


0.5995 

(0.168) 


0.4284 

(0.088) 


0.4512 

(0.180) 


0.4190 

(0.12S) 




90 % 


0.8496 

(0.218) 


0.7911 

(0.241) 


0.5501 

(0.113) 


0.5876 

(0.234) 


0.5517 

(0.174) 


2.0 


80 % 


0.8843 

(0.230) 


0.7974 

(0.240) 


0.5707 

(0.118) 


0.6432 

(0.259) 


0.5944 

(0.290) 




90 % 


1.1030 

(0.286) 


1.0914 

(0.513) 


0.7329 

(0.151) 


0.8376 

(0.337) 


0.7864 

(0.256) 


3.0 


80 % 


1.4552 

(0.589) 


1.5826 

(1.296) 


0.8553 

(0.176) 


1.0345 

(0.419) 


0.9528 

(0.313) 




90 % 


1.8145 

(0.737) 


2.4881 

(1.984) 


1.0983 

(0.227) 


1.3472 

(0.546) 


1.2605 

(0.427) 


4.0 


80 % 


2.3431 

(1.375) 


2.6289 

(2.141) 


1.1399 

(,0.235) 


1.4289 

(0.581) 


1.3138 

(0.437) 




90 % 


2.8401 

(1.447) 


4.3272 

(2.778) 


1.4638 

(0.302) 


1.8610 

(0.757) 


1.7390 

(0.601) 
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TABLE 13 



COVERAGE RATIO (HYPOEXP MODEL) 
N = 50, 0=2/3, pj= 2.0, p 2 = 10.0.p 3 = 2.0 



time C.I (%) 



leneth 



BIN 



NOR 



MLE 



J.K 



B.T 



0.5 80 % 



(.8214) 



90 0/ - 



/0 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



25(.08) 

2491.83) 

26(.09) 

(.T8..87) 

81.03) 
284(.9o) 
81.03) 
(.92. .97) 



25(.08) 
2491.83 
26 .09) 
(.78. .87) 

25(.08) 
267 ,S9) 
81.03) 
(.So, .93) 



0(.00) 

26.0 lj 

29S(.99) 

(. 00 ,. 02 ) 

0(.00) 

4< .0 lj 

296(.99) 

(.00;.03) 



1S6(.62) 
109C.36) 
51.02) 
.31, .42) 



( 



157(.52) 
141(.47) 
21 . 01 ) 
(.41, .53) 



220(.73) 
80(.27) 
01 . 00 ) 
(.22,. 32) 

194(. 65) 
1061.35) 
01 . 00 ) 
1.30..4 1 ) 



1.0 80 % 



(.5788) 



90 0/ - 



/O 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



30(. 10) 
2541.85) 
16(.05) 
(.81..89) 



2SO(.93) 
10 .03) 
(.91. .97) 



30(. 10) 
2441.81) 
26(.09) 
(.77..S6) 

15(.05 

2691.90, 

161.05) 

(.S6..93) 



01 . 00 ) 
179(.60) 
1 2 1 ^. 40 ) 
(.54..65) 

01 . 00 ) 
237(.79) 
63 .21) 
(.74..S4) 



651.22) 
2 19(.73) 
11 (.04) 
(.6S,.7S) 

391.13) 

2521.84) 

91/03) 
(.SO. .88) 



59(.20) 

2251.76) 

1 61.05 ) 
(.71, .81) 

321.1 1) 
261 .87) 
71.02) 
(.83. .91) 



1.5 80 % 



(.3936) 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



191.06) 
2531.84 ) 
28 .09) 
(.S0..S8) 

11 (.04) 
2741.91) 
lo(.0o) 
(.8S..95) 




8(.03) 

264(.88> 

281.09) 

(.84,.92) 

41.01) 

284(.95) 

121.04) 

(.92..97) 



42(. 14 
241 .80 
17 .06) 
(.76..S5) 

17(.06) 
275 .92) 
81.03) 
(.89.-95) 



291.10) 
2401.80) 
31 .10) 
(.75..S5) 

15(.05) 
2651.88) 
201.07) 
(.85. .92) 



2.0 80 % 

(.2652) 

90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



26(.09) 
249 .83) 
25(.0Sj 
(.79,.87) 

51.02) 
283(.94) 
1 21 . 04 ) 
(.92 .97) 



26(.09) 
2331.78) 
41 .14) 
(.73,.82) 

14(.05) 

2611.87) 

25(.08) 

(.83,.91) 



27(.09) 

2611.87) 

12(.04) 

(.83,.91) 

11 (.04) 
285(.95) 
41.01) 
(.93, .97) 




IS(.06) 
274 .91) 
81.03) 
(.88, .95) 



3.0 80 % 



(.1197) 



90 0/ ' 



/ 0 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



141.05) 

27 1(.90) 



(.S7..94 

51.02) 

2921.97) 

3(.01) 
(.96. .99) 



14(.05) 
243 .81) 
43 .14) 
(.77,.85) 

51.02) 

252(.84) 

43(.14) 

(.SO.SSf 



69(.23) 
227 .76) 

41.01) 
(.71, .81) 

376.12) 
261 .87) 

21 . 01 ) 

(.83, .91) 



35(. 12) 
249(.83) 
16 .05) 
(.79, .87) 

17(.06) 
275(.92) 
S(.03V 
(.8b. .95) 



226.07) 
2296.76) 
49(. 16j 
(-72,-81) 

61.02) 
26 16.S7) 
33 .11) 
(.S3..91 ) 



4.0 80 % 



(.0540) 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



136.04) 
275(.92) 
12( .04) 
(.89, .95) 

41.01) 

2961.98) 



(.97,.99) 



136.04) 

2211.74) 

661.22) 
(.69, .79) 

41.01) 
2301.77) 
66(. 22 ) 
(.72, .SI) 



9S6.33) 
1 996-66) 
31.01) 
(.61, .72) 

631.21) 
237(.79) 
06-00Y 
(.74, .84) 



346.1 1) 
248 .83) 
18(.06) 
(. 7S,.S 7) 

151.05) 

2771.92) 

( .m 



196.06) 
229 .76) 
521.17) 
(.72, .81) 

61.02) 
2601.87) 
34(.l lj 
(.83, .91) 
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TABLE 14 



LENGTH OF Cl FOR in PfD> t}(HYPOEXP MODEL) 
N = 50,6 = 2/3^! = 2.0,p 2 = 10.0, p 3 = 2.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.1956 

(0.032) 


0.1704 

(0.031) 


0.1194 

(0.014) 


0.1496 

(0.057) 


0.1351 

(0.035) 




90 % 


0.2450 

(0.040) 


0.2197 

(0.040) 


0.1533 

(0.018) 


0.1983 

(0.076) 


0.1773 

(0.046) 


1.0 


80 % 


0.3463 

(0.049) 


0.3140 

(0.046) 


0.2378 

(0.028) 


0.2267 

(0.069) 


0.2140 

(0.042) 




90 % 


0.4355 

(0.062) 


0.4056 

(0.060) 


0.3054 

(0.036) 


0.3005 

(0.091) 


0.2768 

(0.053) 


1.5 


80 % 


0.5128 

(0.080) 


0.4700 

(0.075) 


0.3563 

(0.042) 


0.3751 

(0.012) 


0.3537 

(0.072) 




90 % 


0.6447 

(0.100) 


0.6117 
(0.100) • 


0.4576 

(0.054). 


0.4972 
. (0.159) 


0.4588 

(0.095) 


2.0 


80 % 


0.7076 

(0.143) 


0.6531 

(0.135) 


0.4748 

(0.057) 


0.5405 

(0.174) 


0.4977 

(0.106) 




90 % 


0.8881 

(0.176) 


0.8630 

(0.193) 


0.6097 

(0.073) 


0.7164 

(0.232) 


0.6584 

(0.144) 


3.0 


80 % 


1.2845 

(0.532) 


1.2267 

(0.526) 


0.7118 

(0.085) 


0.8833 

(0.287) 


0.8240 

(0.182) 




90 % 


1.5970 

(0.591) 


1.8844 

(1.190) 


0.9141 

(0.109) 


1.1708 

(0.380) 


1.0761 

(0.240) 


4.0 


80 % 


2.2724 

(1.102) 


2.6235 

(1.778) 


0.9488 

(0.113) 


1.2309 

(0.402) 


1.1452 

(0.259) 




90 % 


2.7889 

(1.216) 


4.1527 

(2.227) 


1.2184 

(0.145) 


1.6314 

(0.532) 


1.4950 

(0.336) 



39 



TABLE 15 



COVERAGE RATIO (HYPOEXP MODEL) 
N = 20,0-2/3, Pj- 2.0, p 2 = 10.0. p 3 = 2.0 



time C.I (%) 



leneth 



BIN 



NOR 



MLE 



•J.K 



B.T 



0.5 80 % 



(.8214) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



31.01) 
277(.92) 
20(07) 
(. S9..95) 

31.01) 

288(.96) 



261.09) 

2541.85) 

20(.07) 

(.81..88) 

26(.09) 
2651. SS) 
91.03 V 
(.85. .92) 



01.00) 

54(.18) 

2461.82) 
(.14, .22) 

01.00) 

1 15(.38) 
1851.62) 
(.33 .44) 



1241.41) 

172(.57) 

4101) 
(.52. .63) 

881.29) 

211 (. 70 ) 



(.65..76 



1981.66) 

101 (. 34 ) 

11.00) 
(.28, .39) 

144(.48) 
1561.52) 
01.00) 
(.28. .39) 



1.0 SO % 



(.5788) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 




251. OS) 
2421.81 
33(.ll 
(.76, .85) 

221.07) 
2601.87 
1S(.Q6) 
(.83. .91) 



01.00) 

238(.79) 

62(.21) 

(.75,.84) 

0(.00) 
274(.91) 
26(09) 
(.88 .95) 



621.21) 

222(.74) 

161.05) 

(.69,.79) 

371.12) 
2531.84) 
10(03) 
(.80. .88) 



57(. 19) 
2161.72) 
27(09) 
(.67,. 77) 

201.07) 

2611.87) 

191.06) 

(.S3..91) 



1.5 80 % 



(.3936) 



90 0/ - 



/O 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



16(.05) 
257 .86) 
27 .09) 
(.82..90) 

71.02) 

2831.94) 

10(03) 

(.92..97) 



341.11) 
2181.73' 
481.16 
(.68.. 78’ 



161.06) 

2571.86) 

27(09) 

(.S2..90) 



111.04) 
2641.88) 
25(08) 
(.84, .92) 

41.01) 
284(.9 
12(04 
(.92..97' 



2371.79) 
14 .05) 
(.741.84) 

221.07) 
2691.90) 
9(03) 
(.So.. 93) 



25(.08) 
2311.77) 
44 .15) 
(.72.. 82) 

S(.03> 

264(.SS) 

2S(.09) 

(.S4..92) 



2.0 80 % 



(.2652) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



161.05) 
245 .82) 
39 .13) 
(.77,. 86) 

51.02) 

283(.94) 

12(04) 

(.92..97J 



161.05) 

2251.75) 

59(.20) 

(.70..80) 

16(.05) 

2451.82) 

39(.l3) 

(.77..861 



241.08) 
2591.86) 
17(.06) 
(.82,. 90) 

1 11.04) 
2811.94) 
81.03) 
(.91. .96) 



411.13) 

243(.81) 

16(.07) 

(.77,.S5) 

201.07) 
27 1(.90) 
91.03) 
(.87. .94) 



231.08) 
23 11.77) 
461.15) 
(.72,-82) 

71.02) 
2631.S8) 
30(. 10) 
(.33.-92) 



3.0 SO % 



(.1197) 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



201.07) 
257(.S6 
23(.08 
(.82 .90 



31.01) 

2971.99) 

01 . 00 ) 
(.98. .99) 



201.07) 
1891.63) 
91(.30) 
(.58, .68) 

3(01) 

2741.91) 

23(08) 

(.88.95) 



581.19) 
231(.77) 
11(04) 
(.72, .82) 

26(09) 

2691.90) 

51.02) 

(.86,.93) 



391.13) 
242(.81) 
19 .06) 
(.76, .85) 

19(06) 
275 .92) 
6 ( 02 ) 
(.89, .9o) 



18(06) 
2261.75) 
56(.19) 
(.70.. 80) 

61.02) 
262(.87) 
32 .11) 
( .S3..92 ) 



4.0 80 % 



(.0540) 



90 % 



too high 
cover 
too low 
int.val 

too high 
cover 
too low 
int.val 



26(09) 

2741.91) 

OlOOV 

(-8S..95) 

61.02) 

2941.98) 

0 ( 00 ) 

(.96,.99) 



6 ( 02 ) 
1961.65) 
981.33) 
(.60, .71) 

01 . 00 ) 

2021.67) 

9S(.33) 

(.62.73) 



651.22) 
226(.75) 
9103) 
(.71, .80) 

411.14) 

2561.85) 

(.as 



361.12) 
247 .82) 
17(06) 
(.78, .87) 

16(05) 
2781.93) 
6 ( 02 } 
(.90, .96) 



15(05) 
2291.76) 
561.19) 
( .72.. S 1 ) 

61.02) 
2601.87) 
341.11) 
(.S3, .91) 
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TABLE 16 



LENGTH OF C.I FOR In PfD> t)(HYPOEXP MODEL) 
N = 20,0 =2/3,p 1 = 2.0,p 2 = 10.0, p 3 = 2.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.3368 

(0.088) 


0.2706 

(0.088) 


0.1951 

(0.042) 


0.2539 

(0.149) 


0.2382 

(0.105) 




90 % 


0.4186 

(0.110) 


0.3466 

(0.119) 


0.2505 

(0.054) 


0.3306 

(0.194) 


0.3249 

(0.143) 


1.0 


80 % 


0.6008 

(0.158) 


0.5234 

(0.150) 


0.3892 

(0.084) 


0.3918 

(0.161) 


0.3675 

(0.124) 




90 % 


0.7494 

(0.197) 


0.6858 

(0.208) 


0.4998 

(0.108) 


0.5103 

(0.209) 


0.4872 

(0.174) 


1.5 


80 % 


0.9084 

(0.310) 


0.8374 

(0.465) 


0.5840 

(0.129) 


0.6753 

(0.305) 


0.6142 

(0.213) 




90 % 


1.1331 

(0.387) 


1.1701 

(0.S34) 


0.7500 

(0.166) 


0.S794 

(0.397) 


0.8144 

(0.298) 


2.0 


80 % 


1.3743 

(0.845) 


1.3323 

(1.078) 


0.7800 

(0.183) 


0.9871 

(0.470) 


0.8909 

(0.316) 




90 % 


1.6939 

(0.930) 


2.1678 

(2.013) 


1.0016 

(0.235) 


1.2855 

(0.612) 


1.1786 

(0.443) 


3.0 


80 % 


2.5888 

(1.497) 


2.9300 

(2.330) 


1.1776 

(0.357) 


1.6313 

(0.810) 


1.4596 

(0.534) 




90 % 


3. 1 ISO 
(1.550) 


4.7553 

(2.S63) 


1.5122 

(0.462) 


2.1245 

(1.055) 


1.9376 

(0.737) 


4.0 


80 % 


4.2251 

(2.031) 


3.3399 

(2.307) 


1.5894 

(0.763) 


2.2834 

(1.154) 


2.0356 

(0.749) 




90 % 


4.7925 

(1.883) 


4.5249 

(3.430) 


2.0410 

(0.980) 


2.9739 

(1.503) 


2.7125 

(1.03S) 
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IV. CONCLUSIONS 



This thesis considers the problem of estimating the log probability for a semi- 
Markov process which does not enter a particular state before time t. 

Simulation is used to study procedures for obtaining confidence intervals for the 
In probability a semi-Markov process enters a fixed state after time t. The data arise 
from observing N individuals. Three estimators and associated confidence interval 
procedures are compared: The three estimators use different amounts of information 
about the process. The maximum likelihood estimator and its normal confidence 
interval uses the most information. An estimator based on the observed first passage 
times uses the least. An estimator based on an exponential approximation to the 
survivor function of the first passage time uses an intermediate amount of information: 
confidence intervals for this last estimator are obtained using jackknife and bootstrap 
procedures. The simulation results indicate the following. 

(1) Larger numbers of individuals result in smaller average confidence interval 
lengths. 

(2) The binomial confidence interval tends to overcover. This is true in general 
since the target coverage probability is a lower bound on the true coverage 
probability. They also tend to have larger average length. 

(3) The maximum likelihood confidence intervals have the correct coverage if the 
model on which they are based is correct. If the model is incorrect they can 
either overcover or undercover. 

(4) If a jackknife interval does not cover the true value the interval will tend to 
be 'too high'. 

(5) The confidence intervals using the asymptotic renewal estimator have the 
correct coverage for largish t. 

(6) The bootstrap confidence interval requires much more computation than the 
jackknife confidence interval. Since the results for the jackknife and bootstrap 
procedures are similiar, ease of computation appears to favor the jackknife 
procedure. 
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APPENDIX A 

TRUE PROBABILITY TABLE I 

TABLE 17 

EXPONENTIAL MODEL P(D>TIME}. 



: i 

time 


p L = 1.0 .p 2 = 10.0 
6 = 0.5 


Pj = 1-0 .p, = 2.0 
0 = 0.5 


p, = 1.0, p, = 10.0 

0 = 0.6667 


0.0 


1.0000 


1.0000 


1.0000 


0.5 


0.7866 


0.7968 


0.7231 


1.0 


0.6203 


0.6503 


0.5241 


1.5 


0.4891 


0.5351 


0.3799 


2.0 


03S57 


0.4415 


*0.2753 


2.5 


0.3042 


0.3646 


0.1995 


3.0 


0.2399 


0.3012 


0. 1446 


3.5 


0.1S91 


0.2488 


0.1048 


4.0 


0.1492 


0.2055 


0.0760 


4.5 


0.1176 


0.1698 


0.0551 


5.0 


0.0928 


0.1403 


0.0399 


5.5 


0.0731 


0.1159 


0.0289 


6.0 


0.0577 


0.0957 


0.0210 


6.5 


0.0455 


0.0791 


0.0152 


7.0 


0.0359 


0.0653 


0.0110 


7.5 


0.0283 


0.0540 


0.0080 


8.0 


0.0223 


0.0446 


0.0058 


8.5 


0.0176 


0.0368 


0.0042 


9.0 


0.0139 


0.0304 


0.0030 



The sojourn time in state 1 has an exponential distribution with mean 1/pj. The 
sojourn time in state 2 has an exponential distribution with mean 1/p,. 
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APPENDIX B 

TRUE PROBABILITY TABLE II 

TABLE 18 

HYPOEXPONENTIAL MODEL P{D>TIME}. 



time 


p^.O .p : = 10.0 
p^ = 2.0 .9 = 0.5 


Pl = 2.0 .p 2 = 2.0 

p r = 2.0 .9 = 0.5 


p t = 2.0,p : = 10.0 

p, = 2.0.9 = 0.666" 


0.0 


1.0000 


1.0000 


1.0000 


0.5 


0.8652 


0.8670 


0.8214 


1.0 


0.6743 


0.6897 


0.5788 


1.5 


0.5157 


0.5519 


0.3936 


2.0 


0.3930 


0.4465 ' 


0.2652 


2.5 


0.2992 


0.3268 


0.1783 


3.0 


0.2278 


0.2952 


0.1197 


3.5 


0.1734 


0.2402 


0.0804 


4.0 


0.1320 


0.1954 


0.0540 


4.5 


0.1005 


0.1590 


0.0363 


5.0 


0.0765 


0.1293 


0.0244 


5.5 


0.0582 


0.1052 


0.0164 


6.0 


0.0443 


0.0856 


0.0110 


6.5 


0.0337 


0.0697 


0.0074 


7.0 


0.0257 


0.0567 


0.0050 


7.5 


0.0196 


0.0461 


0.0033 


8.0 


0.0149 


0.0375 


0.0022 


8.5 


0.0113 


0.0305 


0.0015 


9.0 


0.0086 


0.0248 


0.0010 



The sojourn time in state 1 is the sum of exponential distribution with mean 1 p i 
and mean l/p.. The sojourn time in state 2 has an exponential distribution with mean 
l,p 2 . 
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APPENDIX C 

AVERAGE LENGTH OF C.I FOR P{D>T} 

Recorded in the Appendix are the average length of confidence intervals for 
P{D> t} corresponding to those given in the thesis for In P{D> t}. If a confidence limit 
is greater than 1, it is set equal to 1. If a confidence limit is less than 0, it is set equal 
to 0. The results are similar to those in the previous tables. 



TABLE IS 

LENGTH OF C.I FOR P[D>t} (EXPONENTIAL MODEL) 
N = 50, 0 = 0.5, p|= 1.0, p 2 = 10.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.1651 

(0.014) 


0.1469 

(0.014) 


0.0692 

(0.007) 


0.1259 

(0.058) 


0.1044 

(0.027) 




90 % 


0.2054 

(0.018) 


0.1887 

(0.019) 


0.0889 

(0.009) 


0.1669 

(0.077) 


0.1330 

(0.033) 


1.0 


80 % 


0.1915 

(0.006) 


0.1741 

(0.007) 


0.1078 

(0.007) 


0.1278 

(0.036) 


0.1154 

(0.015) 




90 % 


0.2389 

(0.008) 


0.2236 

(0.009) 


0.1384 

(0.009) 


0.1694 

(0.048) 


0.1483 

(0.018) 


1.5 


80 % 


0.1967 

(0.002) 


0.1795 

(0.003) 


0.1264 

(0.005) 


0.1381 

(0.035) 


0.1289 

(0.014) 




90 % 


0.2455 

(0.003) 


0.2305 

(0.003) 


0.1623 

(0.006) 


0.1830 

(0.047) 


0.1650 

(0.019) 


2.0 


80 % 


0.1913 

(0.007) 


0.1740 

(0.007) 


0.1320 

(0.002) 


0.1427 

(0.035) 


0.1334 

(0.015) 




90 % 


0.2387 

(0.008) 


0.2234 

(0.009) 


0.1695 

(0.002) 


0.1892 

(0.046) 


0.1711 

(0.019) 


3.0 


80 % 


0.1706 

(0.012) 


0.1526 

(0.013) 


0.1220 

(0.008) 


0.1325 

(0.031) 


0.1233 

(0.016) 




90 % 


0.2123 

(0.016) 


0.1960 

(0.016) 


0.1567 

(0.010) 


0.1756 

(0.041) 


0.1569 

(0.021) 


4.0 


80 % 


0.1459 

(0.017) 


0.1271 

(0.017) 


0.1007 

(0.013) 


0.1092 

(0.028) 


0.1009 

(0.018) 




90 % 


0.1820 

(0.021) 


0.1631 

(0.023) 


0.1293 

(0.017) 


0.1447 

(0.037) 


0.1283 

(0.023) 
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TABLE 19 



LENGTH OF C.I FOR P(D>t) (EXPONENTIAL MODEL) 
N = 20. 0 = 0.5, = 1.0, p,= 10.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.2689 

(0.037) 


0.2262 

(0.045) 


0.1111 

(0.018) 


0.1746 

(0.065) 


0.1584 

(0.055) 




90 % 


0.3281 

(0.046) 


0.2S91 

(0.060) 


0.1426 

(0.023) 


0.2274 

(0.085) 


0.2094 

(0.078) 


1.0 


80 % 


0.3087 

(0.020) 


0.2702 

(0.023) 


0.1708 

(0.018) 


0.1913 

(0.050) 


0.1820 

(0.035) 




90 % 


0.3775 

(0.024) 


0.3469 

(0.030) 


0.2194 

(0.023) 


0.2492 

(0.065) 


0.2356 

(0.047) 


1.5 


80 % 


0.3169 

(0.010) 


0.2792 

(0.011) 


0.1985 

(0.010) 


0.2165 

(0.051) 


0.2030 

(0.029) 




90 % 


0.3876 

(0.013) 


0.3585 

(0.014) 


0.2549 

(0.013) 


0.2820 

(0.-067) 


0.2617 

(0.040) 


2.0 


80 % 


0.3082 

(0.019) 


0.2698 

(0.021) 


0.2058 

(0.007) 


0.2244 

(0.049) 


0.2077 

(0.028) 




90 % 


0.3768 

(0.024) 


0.3463 

(0.027) 


0.2643 

(0.009) 


0.2923 

(0.064) 


0.2650 

(0.037) 


3.0 


80 % 


0.2758 

(0.032) 


0.2345 

(0.037) 


0.1886 

(0.022) 


0.2027 

(0.049) 


0.1854 

(0.034) 




90 % 


0.3365 

(0.040) 


0.3001 

(0.049) 


0.2422 

(0.029) 


0.2640 

(0.063) 


0.2338 

(0.044) 


4.0 


80 % 


0.2377 

(0.045) 


0.1895 

(0.057) 


0.1554 

(0.032) 


0.1627 

(0.050) 


0.1485 

(0.040) 




90 % 


0.2896 

(0.055) 


0.2389 

(0.077) 


0.1992 

(0.042) 


0.21 IS 
(0.065) 


0.1860 

(0.049) 
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TABLE 20 



LENGTH OF C.I FOR P[D> t} (EXPONENTIAL MODEL) 
N = 50, 0 = 2/?. p, = 1.0, p 2 = 10.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.1719 

(0.012) 


0.1601 

(0.012) 


0.0855 

(0.008) 


0.1873 

(0.074) 


0.1645 

(0.044) 




90 % 


0.2160 

(0.015) 


0.2056 

(0.016) 


0.1098 

(0.010) 


0.2482 

(0.098) 


0.2158 

(0.056) 


1.0 


80 % 


0.1962 

(0.003) 


0.1790 

(0.004) 


0.1223 

(0.006) 


0.1523 

(0.042) 


0.1479 

(0.023) 




90 % 


0.2448 

(0.004) 


0.2298 

(0.005) 


0.1571 

(0.007) 


0.2019 

(0.056) 


0.1897 

(0.031) 


1.5 


80 % 


0.1911 

(0.007) 


0.1737 

(0.007) 


0.1320 

(0.002) 


0.1542 

(0.043) 


0.1474 

(0.019) 




90 % 

% 


0.2383 

(0.009) 


0.2231 

(0.009) 


0.1694 

(0.002) 


0.2044 

(0.057) 


0.1907 

(0.025) 


2.0 


80 % 


0.1772 

(0.011) 


0.1595 

(0.012) 


0.1269 

(0.006) 


0.1496 

(0.042) 


0.1397 

(0.019) 




90 % 


0.2208 

(0.015) 


0.2048 

(0.016) 


. 0.1629 
(0.008) 


0.1982 

(0.055) 


0.1803 

(0.024) 


3.0 


80 % 


0.1435 

(0.019) 


0.1245 

(0.020) 


0.0996 

(0.013) 


0.1178 

(0.037) 


0.1073 

(0.022) 




90 % 


0.1780 

(0.024) 


0.1598 

(0.026) 


0.1279 

(0.017) 


0.1561 

(0.049) 


0.1372 

(0.028) 


4.0 


80 % 


0.1140 

(0.021) 


0.0932 

(0.024) 


0.0701 

(0.015) 


0.0812 

(0.031) 


0.0733 

(0.021) 




90 % 


0.1409 

(0.026) 


0.1185 

(0.032) 


0.0900 

(0.020) 


0.1076 

(0.041) 


0.0934 

(0.02S) 
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TABLE 21 



LENGTH OF C.I FOR P(D>t} (EXPONENTIAL MODEL) 
N = 20, 0 = 2/3, pj = 1.0, p 2 = 10.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.2874 

(0.030) 


0.2472 

(0.033) 


0.1374 

(0.021) 


0.3199 

(0.230) 


0.29S4 

(0.158) 




90 % 


0.3508 

(0.037) 


0.3167 

(0.045) 


0.1765 

(0.027) 


0.4166 

(0.300) 


0.4154 

(0.238) 


1.0 


SO % 


0.3151 

(0.012) 


0.2771 

(0.013) 


0.1926 

(0.014) 


0.2595 

(0.098) 


0.2520 

(0.077) 




90 % 


0.3854 

(0.015) 


0.3560 

(0.016) 


0.2473 

(0.018) 


0.3380 

(0.127) 


0.3355 

(0.121) 


1.5 


so % 


0.3057 

(0.022) 


0.2670 

(0.024) 


0.2046 

(0.011) 


0.2539 

(0.076) 


0.2387 

(0.056) 


% 


90 % 


0.3736 

(0.027) 


0.3428 

(0.031) 


0.2627 

(0,014) 


0.3306 

(0.099) 


0.3084 

(0.07S) 


2.0 


80 % 


0.2847 

(0.034) 


0.2440 

(0.038) 


0.2947 

(0.021) 


0.2338 

(0.071) 


0.2140 

(0.051) 




90 % 


0.3475 

(0.042) 


.0.3123 

(0.051) 


0.2499 

(0.028) 


0.3044 

(0.093) 


0.2730 

(0.066) 


3.0 


80 % 


0.2334 

(0.047) 


0.1842 

(0.061) 


0.1515 

(0.036) 


0.1700 

(0.066) 


0.1542 

(0.054) 




90 % 


0.2842 

(0.056) 


0.2316 

(0.081) 


0.1940 

(0.04S) 


0.2213 

(0.086) 


0.1936 

(0.066) 


4.0 


SO % 


0.1893 

(0.051) 


0.1253 

(0.075) 


0.1069 

(0.039) 


0.1117 

(0.057) 


0.1024 

(0.049) 




90 % 


0.2316 

(0.06O) 


0.1536 

(0.096) 


0.1351 

(0.052) 


0.1455 

(0.074) 


0.1281 

(0.060) 
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TABLE 22 



LENGTH OF C.I FOR P(D>t} (HYPOEXP MODEL) 
N- 50,0 = 0.5, pj = 2.0, p2= 10.0, p 3 = 2.0 



time 


C.I <%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.1417 

(0.018) 


0.1217 

(0.019) 


0.0692 

(0.006) 


0.0813 

(0.035) 


0.0722 

(0.018) 




90 % 


0.175S 

(0.022) 


0.1574 

(0.024) 


0.08S8 

(0.008) 


0.1077 

(0.047) 


0.0931 

(0.022) 


1.0 


80 % 


0.1855 

(0.009) 


0.1680 

(0.010) 


0.1078 

(0.007) 


0.1110 

(0.032) 


0.1024 

(0.014) 




90 % 


0.2313 

(0.012) 


0.2157 

(0.013) 


0.1384 

(0.009) 


0.1472 

(0.042) 


0.1317 

(0.019) 


1.5 


80 % 


0.1964 

(0.003) 


0.1792 

(0.003) 


0.1265 

(0.004) 


0.1344 

(0.035) 


0.1245 

(0.016) 




90 % 


0.2451 

(0.004) 


0.2301 

(0.004) 


0.1625 
. (0.006) 


0.1182 

(0.047) 


0.1597 

(0.021). 


2.0 


80 % 


0.1920 

(0.006) 


0.1747 

(0.006) 


0.1322 

(0.002) 


0.1413 

(0.036) 


0.1308 

(0.017) 




90 % 


0.2395 

(0.008) 


0.2243 

(0.008) 


0. 1698 
(0.002) 


0.1872 

(0.048) 


0.1676 

(0.021) 


3.0 


80 % 


0.1675 

(0.014) 


0.1495 

(0.015) 


0.1223 

(0.007) 


0.1265 

(0.032) 


0.1174 

(0.019) 




90 % 


0.2185 

(0.018) 


0.1919 

(0.015) 


0.1573 

(0.009) 


0.1677 

(0.043) 


0.1496 

(0.022) 


4.0 


80 % 


0.1390 

(0.018) 


0.1198 

(0.019) 


0.1009 

(0.011) 


0.0988 

(0.02S) 


0.0919 

(0.019) 




90 % 


0.1723 

(0.023) 


0.1537 

(0.025) 


0.1296 

(0.015) 


0.1310 

(0.037) 


0.116S 

(0.023) 
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TABLE 23 



LENGTH OF C.I FOR P(D> t) (HYPOEXP MODEL) 
N = 20, 0 = 0.5, p j = 2.0, p 2 = 10.0, p 5 = 2.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.2328 

(0.044) 


0.1839 

(0.057) 


0.1108 

(0.016) 


0.1260 

(0.051) 


0.1177 

(0.044) 




90 % 


0.2836 

(0.053) 


0.2314 

(0.077) 


0.1423 

(0.021) 


0.1641 

(0.067) 


0.1592 

(0.065) 


1.0 


80 % 


0.3005 

(0.023) 


0.2614 

(0.026) 


0.1710 

(0.016) 


0.1736 

(0.050) 


0.1633 

(0.032) 




90 % 


0.3671 

(0.029) 


0.3354 

(0.034) 


0.2196 

(0.021) 


0.2261 

(0.066) 


0.2105 

(0.043) 


1.5 


80 % 


0.3163 

(0.011) 


0.2785 

(0.011) 


0.1990 

(0.010) 


0.2110 

(0.053) 


0.1954 

(0.031) 




90 % 


0.3869 

(0.013) 


0.3576 
• (0.015) 


0.2556 

(0.012) 


0.2747 

(0.069). 


0.2502 

(0.040) 


2.0 


80 % 


0.3104 

(0.016) 


0.2722 

(0.017) 


0.2067 

(0.006) 


0.2195 

(0.0501 


0.2017 

(0.032) 




90 % 


0.3796 

(0.020) 


0.3495 

(0.023) 


0.2654 

(0.007) 


0.2859 

(0.065) 


0.2571 

(0.038) 


3.0 


80 % 


0.2709 

(0.035) 


0.2289 

(0.040) 


0.1896 

(0.020) 


0.1917 

(0.047) 


0.1756 

(0.037) 




90 % 


0.3304 

(0.043) 


0.2924 

(0.054) 


0.2434 

10.025) 


0.2497 

(0.062) 


0.2211 

(0.045) 


4.0 


80 % 


0.2276 

(0.046) 


0.1771 

(0.061) 


0.1560 

(0.029) 


0.1470 

(0.047) 


0.1354 

(0.040) 




90 % 


0.2774 

(0.055) 


0.2223 

(0.082) 


0.2001 

(0.038) 


0.1914 

(0.061) 


0.1696 

(0.050) 



50 



TABLE 24 



LENGTH OF C.l FOR P(D> t} (HYPOEXP MODEL) 
N = 50, © = 2/3, p. = 2.0, p^= 10.0, p 3 =2.0 



time 


C.l (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.1560 

(0.015) 


0.1375 

(0.016) 


0.0855 

(0.007) 


0.1349 

(0.050) 


0.1237 

(0.034) 




90 % 


0.1938 

(0.019) 


0.1766 

(0.020) 


0.1098 

(0.009) 


0.1788 

(0.066) 


0.1624 

(0.045) 


1.0 


80 % 


0.1943 

(0.005) 


0.1771 

(0.005) 


0.1226 

(0.005) 


0.1330 

(0.036) 


0.1257 

(0.017) 




90 % 


0.2425 

(0.006) 


0.2274 

(0.006) 


0.1574 

(0.006) 


0.1763 

(0.048) 


0.1615 

(0.021) 


1.5 


80 % 


0.1921 

(0.006) 


0.1748 

(0.007) 


0.1324 

(0.001) 


0.1435 

(0.038) 


0.1346 

(0.014) 




90 % 


0.2397 

(0.008) 


0.2245 

(0.00.8) 


0.1700 

(0.002) 


0.1901 

(0.050) 


0.1722 

(0.019) 


2.0 


80 % 


0.1747 

(0.013) 


0.1569 

(0.014) 


0.1273 

(0.005) 


0.1351 

(0.035) 


0.1257 

(0.015) 




90 % 


0.2176 

(0.017) 


0.2015 

(0.018) 


0.1634 

(0.006) 


0.1791 

(0.046) 


0.1606 

(0.020) 


3.0 


80 % 


0.1327 

(0.019) 


0.1132 

(0.021) 


0.0998 

(0.011) 


0.0957 

(0.028) 


0.0885 

(0.018) 




90 % 


0.1644 

(0.024) 


0.1451 

(0.028) 


0.1281 

(0.014) 


0.1268 

(,0.037) 


0.1123 

(0.022) 


4.0 


80 % 


0.0997 

(0.022) 


0.0769 

(0.026) 


0.0700 

(0.013) 


0.0558 

(0.022) 


0.0544 

(0.016) 




90 % 


0.1230 

(0.027) 


0.0963 

(0.036) 


0.0898 

(0.016) 


0.0779 

(0.029) 


0.0690 

(0.021) 
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TABLE 25 



LENGTH OF Cl FOR ?:D> t} (HYPOEXP MODEL) 
N = 10, 0=2,3. p,= 2.0, P2= 10.0, P 3 = 2.0 



time 


C.I (%) 


BIN 


NOR 


MLE 


J.K 


B.T 


0.5 


80 % 


0.2554 

(0.040) 


0.2112 

(0.047) 


0.1370 

(0.018) 


0.2449 

(0.197) 


0.2360 

(0.138) 




90 % 


0.3113 

(0.049) 


0.2683 

(0.065) 


0.1759 

(0.023) 


0.3189 

(0.257) 


0.3274 

(0.192) 


1.0 


80 % 


0.3126 

(0.015) 


0.2746 

(0.016) 


0.1935 
(0.01 1) 


0.2194 

(0.076) 


0.2105 

(0.056) 




90 % 


0.3823 

(0.019) 


0.3525 

(0.021) 


0.2484 

(0.015) 


0.2857 

(0.098) 


0.2744 

(0.080) 


1.5 


80 % 


0.3084 

(0.020) 


0.2699 

(0.022) 


0.2066 

(0.007) 


0.2298 

(0.060) 


0.2134 

(0.040) 




90 % 


0.3770 

(0.025) 


0.3465 

(0.030) 


0.2653 

(0.009) 


0.2993 

(0.078) 


0.2716 

(0.050) 


2.0 


80 % 


0.2799 

(0.039) 


0.2381 

(0.047) 


0.1971 

(0.017) 


0.2095 

(0.057) 


0.1915 

(0.042) 




90 % 


0.3417 

(0.048) 


0.3044 

(0.062) 


0.2530 

(0.022) 


0.2728 

(0.075) 


0.2403 

(0.051) 


3.0 


80 % 


0.21S0 

(0.047) 


0.1648 

(0.064) 


0.1533 

(0.030) 


0.1414 

(0.052) 


0.1297 

(0.043) 




90 % 


0.2657 

(0.057) 


0.2054 

(0.085) 


0.1965 

(0.040) 


0.1842 

(0.068) 


0.1619 

(0.054) 


4.0 


80 % 


0.1692 

(0.04S) 


0.0964 

(0.074) 


0.1076 

(0.033) 


0.0849 

(0.042) 


0.0794 

(0.037) 




90 % 


0.2080 

(0.056) 


0.1163 

(0.092) 


0.1361 

(0.045) 


0.1105 

(0.055) 


0.0995 

(0.046) 
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